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SUMMARY 


The aerodynamic and mechanical design of a high 
pressure ratio (pr = 8), low mass flow (m = 0,907 kg/s), 
centrifugal compressor, typical of designs that might be 
used for small aircraft and gas turbine engines, has been 
undertaken. The design was subject to the requirements 
that (i) the compressor be of "state-of-the-art" perfor- 
mance with respect to efficiency and range, (ii) the 
inducer tip operate at a relative Mach number of 1.2, and 
(iii) the compressor unit be used as a "workhorse" com- 
pressor together with laser Dopper velocimeter (LDV) 
equipment. 

A Preliminary Design {one-dimensional analysis) 
study examined the trade-offs, for both radial and back- 
ward leaning impellers, with regard to aerodynamic per- 
formance, rotative speed and geometrical dimensions that 
are best suited for use with the LDV equipment and sub- 
ject to the above constraints. State-of-the-art perfor- 
mance dictated designs with backward leaning impellers 
The preliminary design established (i) the rotational 
speed set by the 1.2 inducer tip relative Mach number 
reci'irement and (ii) designs with no splitters so as to 
generate large flow path dimensions relative to the LDV 
measurement resolution volume. 

The final stage configuration designed and fabricated 
is a 75,000 RPM, 19 blade impeller with an axial inducer 
and 30° of backward leaning at the tip. The diffuser is 
a 17 vane, radial divergence, channel diffuser. This 
configuration was calculated to have a 77% total-to- 
static efficiency at the design point with 12% flow range 

from choke-to-surge (range is defined as (1-it /m ■. . ) 

X 100) . surge cnoKe 

Two and/or quasi-three-diraensional aerodynamic and 
stress calculations were performed for the above design. 
Critical speed analyses were performed for the high 
speed rotating impeller assembly. 

The inducer was designed by a method of transference 
from axial compressor technology. The remainder of the 
impeller was designed by using quasi-three-dimensional, 
potential flow techniques iterated with boundary layer 
calculations to control the blade loading distributions 
throughout the impeller. The channel diffuser inlet and 
vane-suction surface were designed using a two-dimensional, 
potential flow/boundary layer calculation procedure. 



A complete mechanical design of the compressor con- 
figuration was undertaken since the compressor unit is to 
be'' tested on a low-speed-of-sound (LSM) compressor test 
facility. Stress and vibration calculations were per- 
formed for the titanium impeller. Adequate burst speed 
was calculated and no problems associated with the blade 
vibration studies were found. However, the finite ele- 
ment disc stress calculations indicate that the bore of 
the impeller will yield locally and stress relieve at 
120% air design speed (90,000 RPM) . 

Apart from the LDV equipment, the compressor has 
been heavily instrumented to obtain the aerodynamic data 
necessary for the analysis of the stage component perform- 
ance. This includes the measurement of total pressures 
in the diffuser and time-averaged and dynamic (unsteady) 
pressure measurements in the inlet, on the impeller cover 
and throughout the vaned, radial diffuser. 

The stage design has been tested for mechanical 
integrity and confirmation of critical speeds up to de- 
sign speed on a companion program. 
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INTRODUCTION 


This report describes the design and fabrication of 
a small mass flow, high pressure ratio centrifugal com- 
pressor stage specifically configured to be used with an 
advanced laser Doppler velocimeter system developed under 
a companion program (NASA/USAAMRDL Contract NAS3-17860) . 
The design of the stage has attempted to be representative 
of the "state-of-the-art” technology and performance of 
the high performance centrifugal used in gas turbine 
engines for helicopter and light aircraft use and off-the- 
road land vehicle applications. At the same time, the 
design has had to meet the practical constraints imposed 
by the use of this hardware with the laser Doppler velo- 
cimeter {LDV) instrument and other diagnostic tools. The 
ultimate purpose of the stage developed on this program 
is its use as a "workhorse" tool for the study of the 
basic fluid dynamics of the inducer, impeller and dif- 
fuser components of the high pressure ratio, centrifugal 
compressor. 

This work is a continuation of the development of 
centrifugal compressors by U. S. government agency spon- 
sored efforts over many years. These development efforts 
date back to the early work by the NACA which was con- 
tinued in the late 1960 's by the U. S. AAMRDL (formerly 
U. S. Army AVLABS) who sparked a surge of centrifugal 
compressor development for advanced helicopter engines. 
Higher efficiency and greater stability to flow distor- 
tions and gas turbine engine transients are the primary 
objectives today for continued research and development 
work on the high pressure ratio centrifugal. Future 
studies using the stage hardware designed and developed 
on this program, in combination with LDV equipment, will 
permit the analysis and improved understanding believed 
to be essential to the development of design tools 
required to achieve the ultimate performance for which 
the centrifugal is believed capable. 

The following sections of this report describe the 
aerodynamic and mechanical design work performed in de- 
veloping this "workhorse" centrifugal stage. 
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DESIGN OBJECTIVES AND SPECIFICATIONS 


The overall objective of this centrifugal com- 
pressor development and fabrication program has been to 
provide a state-of-the-art performance (re to efficiency 
and range) , high pressure ratio, small mass flow centri- 
fugal compressor stage typical of designs that might be 
used for small aircraft and helicopter gas turbine 
engines. 

This stage will be a "workhorse" compressor to be 
used, together v?ith advanced aerodynamic instrumentation 
including, in particular, laser Doppler velocimeter (LDV) 
equipment, on future programs to unravel the complicated 
fluid dynamics within the compressor. The design work 
examined the tradeoffs with regard to aerodynamic perfor- 
mance, rotative speed and geometrical dimensions that are 
best suited for hardware to be used for these aerodynamic 
research purposes. These aerodynamic and physical spec- 
ifications and design constraints are described below, 


Aerodynamic Specifications 

The aerodynamic specifications for a single stage 
compressor (inlet, impeller and matching diffuser and 
collector) have been set down by the contract as follov;s: 


fluid 

inlet ambient temperature 
inlet ambient pressure 

total pressure ratio 
air mass flow 
inlet relative Mach 
number at inducer tip 
no inlet guide vanes 
vane island diffuser 


= air 

= 288. 2°K (518. 7°R) 

= 1.014 X 10^ Newton/ 
meter^ (14.7 psia) 

= 8 

- 0.907 kg/s (2.0 Ib^^^/s) 

= 1,2 

= 0° prewhirl 


Both radial bladed and backswept bladed impellers 
were considered during the preliminary design work. 

A vane island diffuser was to be optimized with 
the impeller design for efficiency and range performance. 

Other Design Constraints 


The ultimate use of the hardware as a "workhorse" 
compressor for aerodynamic studies guided the course of 
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the design. Thus the choices for the final proposed 
design considered the factors of rotational speed and 
md the physical dimensions of the flow passages as they 
impact the use of LDV instrumentation. A low rotational 
speed is desired in order to allow a sufficient time 
period for velocity measurements to be made at a point 
in the flow compared to the time for relative motion of 
the rotating impeller blade passages.* Large dimensions 
of the flow path are advantageous to permit better 
resolution of velocities at a point j.n the flow for a 
given size of the LDV measurement resolution volume. 

The design work included the detailed mechanical 
design of the impeller, vane island diffuser, shroud, 
shaft and other parts required. The mechanical design 
effort included the following requirements: 

1. The impeller should be designed with an 
overspeed Tritrgin of 120% of air design 
speed, to be verified by stress analysis. 

2. The compressor stage mechanical design will 
have no critical speeds between 40% and 
120% of equivalent design speed when tested 
in Creare's Low Speed of Sound Modeling 
facility . 


V: 

This will be primarily achieved in the LDV testing 
by using LSM testing which permits a rotational speed of 
about 0.65 to 0.70 of that required in air testing. 
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PRELIMINARY AERODYNAMIC DESIGN 


The intent of the preliminary design was to set the 
overall geometric quantities and performance levels for 
each component of the compressor subject to the geometric 
and other coiistraints . Both radial and backward leaning 
impellers were analyzed using the standard Creare design 
methodology. Calculations were made using the English 
system of units. 


Design Procedures 


The standard Creare aerodynamic design methodology 
was employed in the investigations of the radial and 
backswept bladed preliminary designs. This design 
methodology is described in detail in Reference 1. 
Additional comments pertaining to specific aspects of 
the preliminary design analysis are discussed below. 

Inlet analysis - An inlet flow analysis and inducer 
optimization study was made to define rotational speed 
and inducer hub to tip radio. This design was conducted 
under the following conditions: (All symbols are defined 

in Appendix A.) 


= 0 (no prewhirl, no IGV s) 

5 2 

Pqq = 1.014 X 10 Newton/meter (14.7 psia) 

Tqq = 288. 2°K (518. 7°R) 

m = 9.072 X 10“^ kg/s (2.0 Ibm/s) 


M 


rel 1* 


1.2 


The following additional specifications were used in 
carrying out the inducer study (these were not sotcified 
under the contract) : 


1} Mass flow is distributed to produce a uniform 
axial velocity profile with radius at the 
inducer inlet (unblocked flow) . 
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2) At design speed, approximately 15% margin in 
mass flow between choke and surge is desired 
(i.e., 7.711 X 10“^ kg/s is the maximum desired 
surge flow rate) . 

3) Blade tip stagger angles of less than 60® are 
desirable at the inducer lead edge. Stagger 
angles larger than 60® make inducer designs 
at high relative tip Mach numbers difficult 
because of excessive blade blockage leading 
to local flow acceleration. 


4) Inducer hub-to-tip ratios 0.3 _< (r|^/r^)^ _< 0.5 

were considered. For each speed N and inducer 
hub-to-tip ratio combination, the 

inducer inlet analysis was used over a range 
of axial inlet velocity C , to determine 

inducer tip relative Mach number and 

inducer tip lead edge relative flow angle, 


A general functional relation N = f(Z^,6j^^) 


exists 


for the given inlet conditions described above. Figure 1 


displays the relation Z 


= (rj^/r_j_) ^ - inducer hub-to-tip 


ratio versus rotational speed N for constant inlet flow 
angle = -65.2° and -60.3° (for it ~ 1-2). 

Figures 2 and 3 display the values of corresponding 

to to the and N range of Figure 1. 


It is clear from Figures 2 and 3 that there are in 
general two solutions (C^j values 1 for which = 1.2. 

For example, at N = 75,000 RPM and = 0.5 these corre- 


spond to 167.6 and 195.7 m/s. The high value of 
corresponds to a low value of The variation of 

at design flow with inlet axial velocity is shown in 

Figure 4 for = 0.5 and in Figure 5 for Z^ = 0.3. 


A second criterion was to maintain a small change in 
flow incidence with mass flow so as to maintain adequate 
inducer stall margin. The change in flow incidence Ai 

m 

with axial velocity at constant inlet relative Mach 

number = 1.2 is shown in Figure 6 for = 0.5 and in 
Figure 7 for = 0.3. 
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Figure 6 " Change in Inducer Tip Incidence 

Between Design Point and Surge 0.5 
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A design axial velocity of 195.7 m/s was selected 

corresponding to an inducer inlet relative flow angle 

3,. = -60.3®. This selection results in a change in 

. ^ . -1 
inlet relative flow angle of 5.25“ between 9.07 x 10 

and 7.711 x 10”^ kg/s (15% flow margin choke to surge). 

A rotational speed of 75,000 RPM at a hub-to-tip ratio 

(r^^/r^)^ = 0.5 was selected to match this relative air 

inlet angle. 


Inducer analysis - Inducer performance has been cal- 
culated using the inducer performance map displayed in 
Figure 8 . This map has been prepared using measured 
transonic axial compressor rotor performance as shown in 
Figure 9 (References 2, 3 and 4). 


The incidence range, choke to stall, at it 

values > 1.0 has been set higher than measured for the 
progenitor axial compressor rotor. Performance measure- 
ments on the axial fan stage indicated that at 90% of 
design speed and above stator performance controlled 
stage stall. On the basis of this evidence, it was 
judged that the axial fan rotor is capable of maintaining 
higher values of unstalled incidence up to M^gi i^ 

than are displayed in Figure 9. This assumption has been 
made in preparing Figure 8 . 


Figure 3 displays the important inducer variables, 
inducer tip diffusion ratio DR as a function of inducer 
tip relative Mach number and inducer tip mean 


incidence angle i^^. To design for a change in incidence 

ancle of Ai = 5.25° between design point flow and surge 
^ m 


flow, the mean incidence i^^^ on the inducer map of 


Figure 8 has been set at 4°. This provides an inducer 
tip diffusion ratio of 1.425 at the design match point 
conditions . 


Impeller design - The Creare design methodology 
utilizes the separated flow model in the impeller and 
incorporates the sudden expansion mixing model in the 
vaneless diffuser space at the exit of the impeller. 

The inducer map and inlet flow analysis described 
above determine the value of diffusion ratio and hence 
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ive Mach number, M^ei it 
Performance Map (lM-70) 






Incidence , ij^(degrees) 


©References 2 and 3 
^Reference 4 






impeller tip jet relative Mach number and relative veloc- 
ity. Several other parameters must be selected for the 
preliminary design analysis. 

The mass flow fraction in the wake has been set 

at 0.2 for all cases analyzed. Experimental measurements 
{e.g., Reference 5) have verified, at least for radial 
bladed impellers, that the mass flow fraction in the wake 
can be expected to vary between 0.15 and 0.25. 

The swirl parameter (= tangent of the swirl 

angle measured to radial at impeller exit) has been set 
at 2.75. It is believed that for high pressure ratio 
centrifugal stages, this is close to the optimum value. 


Backflow from the diffuser vaneless space into the 
impeller is often a serious design ill for high pressure 
ratio compressors. Unfortunately little is known defini- 
tively about this phenomena. 

Our present analytical capabilities are insufficient 
to quantify the expected magnitude of this loss as a 
function of compressor design variables (X 2 */ number of 

diffuser channels, etc.). Throughout the preliminary 
design task a constant backflow loss of three points in 
stage efficiency has been assumed. This may be treated 
as "insurance" on efficiency prediction since recent 
experience with impellers and diffusers designed by Cieare 
methods indicates zero backflow loss. We are thus hope- 
ful that up to three points gain in efficiency above the 
predictions herein might be achieved. However, we are 
reluctant to credit this margin at this time and prefer 
to treat it as reserve and hedge against uncertainties in 
the design predictions. 

Channel diffuser design - Clearly defined single 
plane divergence diffuser performance maps are available 
only for channel diffuser aspect ratios of 0.25, 1,0, and 
5.0 (References 6 and 7). To decide what the optimum 
geometry is for a diffuser aspect ratio not equal to those 
for which performance maps are available, some interpola- 
tion scheme is needed. Such a scheme has been used in 
this preliminary analysis and is described below. 
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The diffuser throat aspect ratio (AS^ = ^ 4 /^ 4 ^ 

been set by requiring b^ “ ^ 2 ' resulting aspect 

ratio for all designs studied lies between 0.25 and 1.0. 
Since there are no geometrical limitations on the value 
of the diffuser exit radius, radial plane (rather than 
meridional plane) diffusion has been employed for the 
channel diffusers. At the high match point blockages 
encountered in this design, there is a distinct possibil- 
ity of achieving better channel diffuser recovery with 
radial plane diffusion rather than meridional plane 
diffusion. The reason stems from the fact that diffuser 
recovery falls off faster for throat aspect ratios 
greater than 1.0 than for throat aspect ratios "less than 
1.0 at high blockages. The situation is reversed for 
low blockages. 

Channel diffuser recovery charts are needed 

for aspect ratios between 0.25 and 1.0. These were 
created for AS^ - 0.25 and 1.0 for a fixed diffuser geom- 
etry of 20 = 10° and L/W = 17. This geometry has been 
chosen as a compromise made between (i) providing a 
sufficient buffer region between this geometry and the 
steep fall off in Cp that occurs at high values of 

20 and (ii) on the other hand not providing too large a 
buffer region and unnecessarily penalizing channel dif- 
fuser recovery and, therefore, stage performance. 

The number of diffuser channels has been set at 17 
to achieve adequate flow area at the channel diffuser 
throat yet maintain AS. > 0.25 for radial plane diver- 
gence , 

Blade stress calculations - There are several 
uncertainties in the stress analysis of blading for cen- 
trifugal compressors. Blading stresses arise from at 
least four different sources; 

1 ) centrifugal loads, 

2 ) vibration, 

3) displacement of boundaries, i.e. , distortions 
at the blade attachment point to the disc due 
to combined thermal and centrifugal stresses 
in the disc, and 

4) thermal stresses within the blade. 
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To calculate all these stresses is impractical; the 
number of boundary conditions is formidable, and a (prop- 
erly manufactured) blade usually fails by frtigue which 
means a ciimulative calculation of the stress history of 
the material is needed. The procedure adopted by most 
stress analysts is to model the loadings on the blade in 
a very simple way. These models are, of course, aug- 
mented as often as possible by insight gained from hard- 
ware tests and/or numerical solutions for simplified 
boundary conditions, 

A fairly standard approach is outlined below for the 
preliminary estimate of blade stresses for radial blades. 
Some special calculations performed for backward leaning 
blading are described later. The approach used for 
radial blades was as follows: 


1) The blades were considered to be cut into 
radial strips to calculate the average 
tensile stress in the root. This calcu- 
lation was performed ignoring the fillet 
radius at the blade root. 


The material properties were derated to 
account for vibration. The diminution of 
properties depends on the type and number 
of blade natural frequencies present with- 
in the operating range of the machine. 
Typical vibration allowances are between 
1.379 X 1 q 4 and 5.516 x 10^ Newton/meter^ . 


3) The impeller disc was designed so as to 
minimize the blading deformations. That 
is, no direct account of the stresses due 
to displaced boundary conditions was 
attempted . 

4) Thermal stresses were estimated at stations 
along the blade from a one-dimensional heat 
transfer analysis. 

The calculations were performed in full only for 
several very critical stress configurations. Usually 
then, only Item (1) was effected for a design and this 
quantity was calculated only because it is necessary 
input to calculate disc stresses. 
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Creare has no prior experience calculating stresses 
in backward leaning impellers/ so a preliminary investi- 
gation of the effects of the backward lean per se was 
undertaken. For this purpose, the centrifugal force was 
decomposed into components along and perpendicular to the 
b].ade. To first order, the increase in stress due to the 
backward lean is contained in the perpendicular components 
of force.* The stresses in the blade were then calculated 
as follows: 

A section was cut perpendicular to the blade and 
extending into the depth of the blade axially. The forces 
due to the centrifugal loading on a small element were 
determined by the mass of the element, hence the density 
of the material, the square of the rotation speed, etc. 
These forces produce moments in the blade section. These 
resultant stresses were calculated using the slender beam 
approximation : 



where ; 

M = moment at the section 

— = section modulus 
c 

It is then a straightforward process to determine the 
distribution of stress along the blade. This has been 
done for two simple shapes represented by a linear and a 
parabolic taper. The variation in stress along the beam 
is different for the two cases; for present purposes, it 
was of interest to compare the maximum stress occurring 
in the root for the two shapes versus taper ratio — 
taper ratio being the ratio of thickness at the root 
compared to the thickness at the tip of the blade. The 
effect of a blade root fillet was not considered in these 
calculations in the anticipation that the calculations 
are conservative. 


*The uniaxial equivalent stress of the biaxial 
stress field is approximately the difference of the 
principal stresses. Thus in the worst case, the com- 
pressive stress due to bending of the blade would add 
to the tensile stress to determine the equivalent stress. 


21 



Preliminary Design Analyses 

The Creare design methodology was applied to a num- 
ber of radial bladed and backswept bladed designs 
(0° 1.^2b — Various inducer performance maps and 

diffuser designs using both meridional and radial plane 
channel diffuser divergence were studied. Only the final 
design calculations using the "best" projected inducer 
performance and radial plane diffuser divergence, as dis- 
cussed in the previous section, are presented here. 

Speed and inducer hub-to-tip ratio - The future use 
of this stage as a "workhorse" compressor dictated the 
choice of a low rotational speed. An upper limit of 
inducer lead edge selected. As des- 

cribed in the previous section, the corresponding rota- 
tional speed of 75,000 RPM was then specified to meet 
the requirement of inducer relative Mach number 
^rel It ~ 1*2. The specific speed of the compressor 

( 82 b “ -30°) at this rotational speed is = 0.633. 

It is common belief that best efficiency performance 
for high pressure ratio compressors is attained at high 
specific speeds (order of = 0.775. The preliminary 

design analysis was used to examine the predicted stage 
efficiency increase at the maximum allowable rig speed 
N = 82,650 RPM for = -65.2° and -60.3°, 

and ~ 0-3. The results for both values 

(corresponding to the values) are tabulated below. 

No great significance should be attached to the differ- 
ence in efficiencies for N = 75,000 RPM cases. 


TABLE I - EFFICIENCY VARIATION WITH SPEED 


Inlet Axial 
Velocity 

(meter s/s) 

Total to Static Efficiency, rt 
(Backward Leaning Angle 82 b ~ -30°) 

Rotational Speed 
N=75,000 RPM;2 o=0.5 

Rotational Speed 
N=82,650 RPM;2 o=0.; 

167.6 

76.3% 

76.8% 

195.7 

76.4% 

76.8% 
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The increase in n of approximately 0.5 points 

T” S 

comes principally from a reduction in friction work. 

This stems from a smaller impeller tip radius (r 2 ) for 

the higher speed case (N =0.698) in providing an equal 

pressure rise across the impeller. The sensitivity of 

efficiency at the design point (n ) is 9 n.,^/9 N = 0 . 6 % 

DP DP 

per 10,000 RPM. 

Thus going to a higher rotational speed (N = 82,650 
RPM) would raise the predicted stage efficiency about 
0.5 points? this was not considered warranted relative 
to the desire to maintain low rotational speed. 

Design point inducer performance/impeller exit 
swirl parameter - Early in the preliminary analysis, a 
study was made of the sensitivity of stage performance 
(Oip^g) and flow path dimensions to the independent design 

parameters of inducer tip diffusion ratio DR and impeller 
exit swirl parameter ^ 2 *- This study was carried out 

over the range of backswept blade angles 0 ^ 62)3 S “40°. 

The results are presented in Figures 10 through 12. 

Figure 10 displays the stage efficiency nip _5 at the 

design point a” a function of inducer DR for &2'b ~ “30° 

and -40° for several values of ~ 2.25 and 2.75). 

The influence of DR is readily apparent as is the effect 
of increased blade backsweep in raising stage efficiency. 

Sensitivity factors are displayed in the table below 

for the match point efficiency with respect to the input 

parameters 6 _, , DR, (m /m) , X^., and M. These sensitiv- 
2 b w 2 * 

ity factors are useful for deciding the optimum set of 
trade-offs in a situation requiring a compromise among 
the several input parameters. 

The sensitivity factors are expressed as 3 Oj^p/9 X 

where X is an independent (input) variable and n^p is the 

design point total-to-static stage efficiency. These 
sensitivity factors were obtained from the results of a 
large number of calculations of design configurations 
using several values of the above input parameters. To 
obtain these factors, two large numbers were subtracted 
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Diffusion ratio, DR 


Figure 10 — Voriotion of Compressor Efficiency with Design Parameters 






from one another and this introduces some uncertainty. 
Further error may result just from the tolerances allowed 
within the computer programs to match the required input 
quantities . 

The following table displays a representative value 
of these sensitivity factors that are valid around the 
range of designs reported herein. 


TABLE II - SENSITIVITY FACTORS 


Rate of 

Change of 

Total-to- 

Static 

Efficiency 

per Change 

in 

Variable , X 
9 n^p/S X 


Variable X 


Backward 
Leaning 
Angle , 

3 

Diffusion 
Ratio f 
DR 

Wake Mass 

Fraction , 

m /m 
w 

Swirl 

Parameter, 

^2* 

Speed, 

N 

0.6 pts 
per -10° 

1 . 6 pts 
per 0.1 

1 

-1.3 pts 
per 0 . 1 

-1.1 pts 
per 0.5 

0.6 pts 
per 
10,000 
RPM 


Figure 11 displays the impeller exit tip depth i >2 as 
a function of DR and I 2 * with ^2h parameter. The 

strong dependence of impeller tip depth on blade angle 
02 J 3 and A 2 * apparent. 

Figure 12 shows the diffuser vaneless space radial 
dimension as a function of the same variables used in 
Figure 11. The vaneless space radial dimension (Ar) 
equals - r 2 , where r 3 = 1.05 x r 2 with r^ defining the 

diffuser vane leading edge radius. The Ar variation is 

primarily a function of blade angle 6_, . 

2b 

These results were used to assist in the selection 
of the proposed ^ 2 * amount of blade backsweep. 

Because of the strong sensitivity of b 2 with change in 

^ 2 *’ ^ value of ~ 2.75 was selected for all of the 

final design calculations . This provided the largest 
value of b 2 consistent with empirical design correlations 
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Vaneless gop, AT = r 3 -r 2 (cms.) Voneless gap.A^” =*3 — ^(cms.) 
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Swirl parameter, X^* = (Cq 

Figure 12 — Voneless Space Radial Dimension Variotion with Design Parameters 

(m^^/m =0-2 , T3= 1.05 r^) 
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on the optimum value of ^ 2 * discussed under "Design Pro- 
cedures". A large value of b 2 is advantageous for use with 
LDV instrumentation. 

DR values have very little influence on the flow 
path physical dimensions. Hence DR at the design point 
was selected to provide reasonable range to choke at 
design speed using the final inducer map selected for 
this design study. The inducer map and selection of 
design DR were discussed previously. 

It is evident that efficiency can be increased with 
backsweep. However, the impeller tip depth b^ dimension 

size is decreased somewhat with large backsweep 
f^2b = -30“ to -40“ ) and, as discussed previously, the 

design proposed was one of some degree of backsweep in 
order to seek increased efficiency performance. 

Disc and impeller blade stress calculations - Using 
the preliminary stress analysis procedures for blade 
stress presented previously, the following calculations 
were made for disc and blade stresses. 

Past experience with many designs has shown that it 
is possible (with careful design) to limit peak disc 
stresses to 

%ax = <“■” ' ^ 


where : 

p = material density 

U 2 = tip speed 

V = Poisson’s ratio 

A preliminary decision was made to use titanium 
6A1-4V for the design because it offered the best 
strength to weight ratio of candidate materials whose 
strength to weight ratios versus temperature were known. 

The table below shows the expected tip speed limits 
for the impeller at a maximum bore temperature of 93 “C. 
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TABLE III - TITANIUM STRESS LIMITS 

Peak Allowable Tip Speeds 

Using 6A1- 

for the Impeller at 93.33°C 
-4V Titanium 

Stress Limit Condition 

Range of Limiting Tip Speeds 
m/s 

0.2% yield stress - 
6.619 X 10° Newton/m^ 

603.5 - 670.6 

Ultimate strength = _ 

7.446 X 10® Newton/m'^ 

640.1 - 746.8 


In suininary, it is feasible to operate at tip speeds 
between 600 and 670 m/s, and with careful design, it 
should be possible to achieve tip speeds as high as 
750 m/s — with well-known materials. It was concluded 
that no problem should arise in designing an impeller to 
run at the highest tip speed proposed to meet the aero- 
dynamic performance targets, provided blade backsweep 
was not too large. 

The following table shows the comparison between the 
non-dimensional stress in the root of the blade and the 
taper ratio for two cases, a linear and a parabolic taper. 


TABLE IV - STRESS. FUNCTION VS. TAPER RATIO CONSTANT 


k 

linear 

0 

0.5 

1.0 

2.0 

3.0 

4.0 

6.0 

taper 
f 2 , para- 
bolic 

0.5 

0.259 

0.167 

0.093 

0.063 

0.047 

0.031 

taper 

0.5 

0.241 

0.146 

0.074 

0.047 

0.033 

0.020 


where : 

t = 1 + kx (linear taper) or 
t = 1 + kx^ (parabolic taper) and 

t = blade normal thickness distribution (normalized) 

X = distance along blade from tip to root (normalized) 
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The quantity f is defined by 

X 

f (x) = [l/t(x)]^ / (x-z)zdz 

o 

and k is related to the taper ratio by k = (t. . /t^ . )-l. 

■* hub^ tip 

Assuming a taper ratio of 2, a normalized stress of 
approximately 0.15, the peak stress in the blade root was 
calculated. The equation for peak stress is: 


fh^ 6pr(jj^ sinS]3 



The preliminary aerodynamic design only gave the blade 
height at the tip of the wheel. The stress was calcu- 
lated below for this case: 

r = ^2 ^ 7.976 x 10“2 m 

N = 75,000 RPM 

Sfa - -30” 

h = ^ 3.302 X 10"^ m 

t , . = 1.27 X 10“® m 

tip 

k = 1.0 

a = 8.481 X 10^ Newton/meter^ 

The above stress of 8.481 x 10^ Newton/m^ was low 
enough that the extra loading due to the centrifugal 
effect on backward curved blades did not appear to be of 
paramount importance. However, it was expected that the 
maximum bending stress would occur at less than the tip 
radius . 

Results 

Aerodynamic states and stage geometry - Table V 
displays the important aerodynamic flow states and stage 
geometry at the design point for four values of impeller 
blade backsweep (32b ~ -20°, -30°, -40°). 

Figure 13 displays the increase in predicted stage 
design point efficiency ri^^p as a function of blade angle 

®2b- 
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TABLE V - RESULTS FOR RADIAL AND BACKSWEPT BLADE 
DESIGNS (DESIGN POINT CONDITIONS) 



Radial Blade 
Design 

Backswept Blade 

Design 

^2b 

0° 

-20« 

-30“ 

-40“ 

Pr 

8 

8 

8 

8 

m 

9.072x10"^ 

9.072x10"^ 

9.072x10"^ 

9.072x10 

N 

75,000 

75,000' 

75,000 

75,000 

M 1 - 10 . 
rel It 

1.2 

1.2 

1.2 

1.2 


195.7 

195.7 

195.7 

195.7 

^It 

343.2 

343.2 

343.2 

343.2 


0 

0 

0 

0 

^It 

-60.3 

-60.3 

-60.3 

-60.3 

PQO 

1.014x10^ 

1.014x10^ 

1.014x10^ 

1.014x10 

T 

oo 

288.2 

288.2 

288.2 

288.2 

^Itb 

-56.3 

-56.3 

-56.3 

-56.3 

^It 

0.583 

0.583 

0.583 

0.583 

^It 

395.3 

395.3 

395.3 

395.3 

^2J 

577.0 

525.8 

498.0 

463.9 

c 

^T2J 

489.8 

445.6 

425.5 

403.6 

*^M2J 

304.8 

279.3 

258.9 

229.2 

c 

^2W 

589.5 

586.1 

585.2 

578.8 

^T2W 

582.5 

580.3 

579.4 

573,0 

*'M2W 

90.2 

82.6 

81.4 

82.75 

^2J 

318.5 

324.3 

327.7 

331.0 


Continued. . 
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Table V - Continued 


^2b 

0° 

-20° 

-30° 

-40° 

«2w 

90.2 

87.87 

94.00 

108.1 

“2J 

1.427 

1.280 

1.201 

1.109 

«2w 

1.356 

1.319 

1.300 

1.271 

^rel 2j 

0.791 

0.791 

0.792 

0.792 

W 

sep 

277.3 

277.3 

277.3 

277.3 

DR 

1.425 

1.425 

1.425 

1.425 

e 

0.494 

0.497 

0.483 

0.450 

^02 

1.117x10^ 

1.034x10® 

9.997x10^ 

9.515x10^ 

p 

o2W 

1.007x10^ 

1.082x10® 

1.138x10® 

1.172x10® 

^02 

1.096x10^ 

1.041x10® 

1.027x10® 

9.997x10^ 

^2 

3.385x10^ 

3.833x10^ 

4.123x10^ 

4.433x10^ 

^2j 

406.2 

420.8 

429.4 

438.3 

*^2W 

470.3 

492.4 

516.7 

518.1 

^o2j 

568.3 

555.3 

550.1 

542.9 

'’•’o2W 

637.2 

656.7 

668.9 

677.8 

^2 

582.5 

610.5 

626.4 

642.2 

^2* 

440.3 

453.4 

460 

465.4 

^o2* 

582.2 

575.6 

573.9 

570 

^o2** 

593.9 

588.3 

587.2 

583.9 

^2* 

3.827x10^ 

4.220x10^ 

4.447x10^ 

4.661x10^ 

^o2* 

1.034x10® 

9.928x10^ 

9.860x10^ 

9.653.10^ 

^2* 

1.288 

1 • X8X 

1.133 

1.080 

^Po2 

5.716x10^^ 

4.806x10'^ 

4.033x10^ 

3.034x10^ 


Continued. . 
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Table V - Continued 


^2b 

O'* 

-20° 

-30° 

-40° 


508. 

472.4 

456.3 

437.4 

^M2* 

184.9 

171.8 

166.0 

159.1 

C2* 

540.7 

502.9 

485.5 

465.4 

^2* 

2.75 

2.75 

2.75 

2.75 

a 

0.873 

0.877 

0.882 

0.889 

M4 

0.95 

0.95 

0.95 

0.95 

^4 

5.709x10^ 

5.571x10® 

5.523x10® 

5.426x10® 

^o4 

1.020x10® 

9.860x10® 

9.791x10® 

9.653x10® 

®4 

0.160 

0.135 

0.118 

0.103 

^2*-4 

0.288 

0.236 

0.199 

0.153 

*"p4-c 

0.541 

0.581 

0.610 

0.634 

'^pD 

0.657 

0.675 

0.685 

0.689 

Angp(%) 

3.0 

3.0 

3,0 

3.0 


0.99 

1.33 

1.56 

1.82 

Aripc (%) 

1.19 

1,60 

1,87 

2.18 

Arigp {%) 

2.19 

2.93 

3.42 

4.00 


2.76 

2.44 

2.08 

1.61 

Arip (%) 

12.04 

10.36 

9.52 

8.87 

Ari3-^(%) 

4.81 

5.00 

5.01 

4.91 

’^ITT 

87.24 

86.62 

86.49 

86.48 

n,pg (%) 

75.20 

76.27 

76,98 

77.62 


Continued. , 
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Table V - Concluded 


^2b 

0° 

Geometry 

-20® 

O 

o 

1 

J 

O 

o 

I — i 

u 

2.187x10“^ 

2.187x10"^ 

2.187x10"^ 

2.187x10“^ 

^It 

4.371x10"^ 

4.371x10"^ 

4.371x10"^ 

4.371x10"^ 

^2 

7.417x10"^ 

7.772x10“^ 

7.976x10"^ 

8.179x10"^ 

b2 

3.480x10"^ 

3.327x10“^ 

3.226x10"^ 

3.175x10"^ 

^4 

3.480x10"^ 

3.327x10“^ 

3.226x10“^ 

3.175x10"^ 

^4 

1.090x10"^ 

1.133x10“^ 

1.148x10"^ 

1.166x10"^ 


0.32 

0.29 

0.28 

0.27 


7.798x10"^ 

8.153x10"^ 

8.382x10"^ 

8.585x10"^ 

Ar(Er3-r2 

} 3.708x10“^ 

3.886x10”^ 

3,988x10"^ 

4.089x10"^ 

#Dif fu- 

sers 

17 

17 

17 

17 

# Blades 

19 

19 

19 

19 


Notes : 




1. 

Mass flow in kilograms per second. 

2. 

Velocities in meters 

per 

second. 

3. 

Pressures in Newtons 

per 

square meter. 

4. 

Temperarures in °K. 



5. 

Angles in degrees . 



6. 

Lengths in meters. 
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Figures 14 and 15 display the predicted performance at 
100% speed for the designs with -30®, respectively. 

Figures 16 and 17 show the inducer inlet and impeller 
exit velocity diagrams, respectively, for the two designs 
displayed in Figures 14 and 15 (02b ” -30°). 

A definition of the impeller meridional plan form 
and detailed blade shapes was not available as this type 
of calculation was performed as part of the final aero- 
dynamic and mechanical design. An estimate of the 
impeller meridional configuration based on prior "experi- 
ence" and assumed blade thickness distribution was used 
in performing the preliminary stress analyses. 

Preliminary stress analysis - There are four chief 
areas of concern regarding the mechanical design of an 
impeller for high tip speed. 

1) The disc material and design must be adequate 
to withstand the centrifugal loading imposed 
upon it during its required lifetime. 

2) The blades must be adequate to withstand the 
centrifugal and thermal loading. 

3) The blades (and disc) should be designed so 
that there are no severe vibration modes which 
can cause failure. 

4) The impeller deflections during operation 
should be smaller than the clearances. 

The preliminary stress analysis was directed toward 
areas (1) and (2) of the above. The aerodynamic design 
had not progressed to the point where it was possible to 
design mechanically with respect to areas (3) and (4) . 

The disc stress analysis indicates that it will be 
possible to operate the proposed designs up to tip speeds 
of order 670 m/s using titanium (6A1-4V) . With careful 
design, tip speeds up to order of 750 m/s should be 
possible. 

The first order blade stress analyses (see section 
on disc and impeller blade stress calculations) indicated 
that the backward leaning blade designs proposed would 
not meet difficulties re to blade stresses. However, 
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these were very preliminary and approximate calculations 
and the final value of 3213 used for the final design had 

to await a more complete blade analysis that was made 
during the final mechanical design. At the same time, 
empirical evidence has demonstrated that impeller designs 
in titanium have been routinely run with 32 jj ~ -32“ at 

tip speeds up to 701 m/s. Therefore, the proposed design 
for 32 J 3 “ -30® appeared to be reasonable of attainment 

at the 120 % overspeed required. 

Conclusions/Recommendations 


The following conclusions/recommendations were made 
on the basis of the preliminary design analysis and the 
contract requirements. 

1. A rotational speed of 75,000 RPM with an 

inducer inlet hub-to-tip ratio = 0.5 was 

proposed. This speed produces a stage 

specific speed N = 0.633. The highest speed 

s 

deemed practical* for this design (N = 82,646 
RPM at hub-tip ratio = 0.3) would produce a 
specific speed = 0.698, but only a predicted 

increase in stage efficiency of order of 0.5 
points. 

2. Stage efficiency increases on the order of 2 
to 2,5 points at the optimum impeller design 
is changed from radial (3233 = 0 °) to backswept 

blades {^2h ~ -40°) respectively. The 

predicted efficiency of a radial bladed impeller 
(n^g = 75.2%) was judged too low to be repre- 
sentative of state-of-the-art performance at 
the design pressure ratio and mass flow. The 
performance (n,j,g = 77.0% for 62 ^^ “ -30® and 

riipg ~ 77.6% for 8233 = -40°) with backsweep 

was sufficiently high to recommend that a back- 
sweep design of 62 J 3 = -30° be selected for the 

final design. 


This upper speed limit is set by the LSM facility 
capability. 
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3 . 


A 19 blade inducer/ impeller was recommended with 
an inducer configuration established on the 
basis of empirical transonic axial compressor 
aerodynamic performance, A design with no 
splitters will reduce the uncertainties of aero- 
dynamic design and provide a design more com- 
patible with the use of LDV equipment. 

4. A 17 channel, radial plane divergence diffuser 
was proposed to optimize the overall diffuser 
pressure recovery Cpp and provide the predicted 

stage efficiencies. 

5. At the impeller tip speed for the proposed back- 
swept design {of -30°), hub and blade stresses 
were not judged to be critical limiting factors 
at 100% speed. However, at 120% overspeed, 
these calculations showed that hub stress may 

be critical for the proposed backswept design. 
Confirmation of the preliminary stress analysis 
was made after the final mechanical design 
analysis. 


42 



AERODYNAMIC DESIGN OF INDUCER 


General Philosophy 

The general impeller design approach followed here 
recognizes two parts of the rotor - an inducer section 
and the remainder of the impeller. This view is being 
updated continuously at Creare as more and more is 
learned about the flow physics of the centrifugal rotor. 
Presently, it appears that the impeller (i.e., complete 
rotor) flow is best divided into an inlet flow (the 
inducer) and an internal flow within the impeller. The 
inlet flow is controlled by incidence and Mach number 
effects. The performance of the inlet region of the 
impeller is strongly determined by blade shape, blade 
thickness, blade number, taper ratio, hub-tip ratio, etc. 
This inlet region may be considered to extend about one 
throat width downstream of the throat. 

Farther into the impeller, the flow appears to 
behave as if it were in a complicated diffuser. Overall 
diffusion within the impeller is controlled by the 
inducer inlet (as mentioned above) , and impeller para- 
meters such as blade shape, blade number, blade loading, 
discharge flow angle, etc. These parameters determine 
the relative velocity distributions on the four surfaces 
bounding the impeller channel. These relative velocity 
distributions and the incidence and relative Mach number 
effects at the inducer inlet determine the capability of 
the wall boundary layers to negotiate the adverse pres- 
sure gradients of the diffuser-like flow in the impeller. 
The boundary layer behavior is also exacerbated by blade 
tip leakage and rotational effects (Coriolis effects) 
which are not well understood. It appears from available 
empirical data that there may be a limiting impeller 
relative velocity diffusion ratio for a given area ratio 
impeller. Just how the impeller aerodynamic and geo- 
metric parameters will control the diffusion ratio with- 
in the impeller is not well established today either by 
theory or experimental evidence. 

As will be described below, the inlet part (i.e., 
the "inducer") of the impeller was designed by a method 
of transference from axial compressors. The centrifugal 
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impeller will then have ensured excellent and predictable 
inlet performance. 

The rear portion of the impeller has been designed 
to give good "diffuser-like" characteristics by con- 
trolling the velocity distributions on critical surfaces 
(the suction surfaces at the tip and hub). By this means, 
the design should attain the impeller discharge conditions 
predicted in the Preliminary Design Phase. 

Inducer Design Method 

The design of the inducer blade geometry is based 
upon a controlled scaling of a high performance transonic 
rotor. The progenitor axial fan used in this design is a 
Pratt and Whitney propulsion fan (Reference 2) . The 
inducer map u'^ed in the Preliminary Design Phase was 
based on data for this transonic rotor. The basis of 
the design procedure is the logical assumption that the 
inducer performance will be replicated if the geometrical 
and aerodynamic parameters of the original transonic 
fan blading are duplicated. Clearly, not all these 
parameters can be duplicated; nevertheless, an 
attempt was made to duplicate what are believed to be the 
most important ones. 

Design parameters - The most important parameters 
are deemed to be : 

- approach relative tip Mach number, 

- blade tip pitch to chord ratio, (s/c)^, 

- approach incidence, i (camber line) , 

m 

- blade suction surface profile, 

- blade "throat" to approach flow area ratio, 

A . /A , . 
mxn rel 

The blade-tip flow region is transonic and is 
believed to be the most critical aero dynamic ally (in 
terms of the above parameters) for duplicating the fan 
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performance. Consequently, emphasis was placed on the tip 
section in performing the basic scaling. The inducer 
blade shape, particularly in the hub region, will diverge 
significantly, behind the throat, from the original trans- 
onic blading of the progenitor fan. This will be evident 
from the design procedures described below. However, as 
a consequence of the low approach relative Mach numbers, 
this portion of the blading is thought to be less critical 
to the replication of the rotor performance than is the tip 
section. 

Finally, the rotor performance is assumed to be 
describable in terms of blade element performance. The 
rotor performance should be essentially duplicated if the 
critical parameters (listed above) of each blade element 
are maintained, even if the blade elements are stacked 
differently than for the original fan. 

The above philosophy is reflected in the detailed 
design procedures below. The two major steps in this 
procedure provide : 

1) fan parameters (i , A . /A , , suction surface 

m min rel 

profiles) from the progenitor blading, and 

2) the blade characteristics for the present 
impeller by the maintenance of these parameter 
values . 

The geometric quantities discussed above are defined in 
Figure 18. 


Design Procedures 

The fan blade profiles are defined in coordinate 
form on a number of cylindrical sections. The relevant 
parameters for each of these sections are to be calcu- 
lated and the constancy of these parameters maintained 
in defining the blading for the present design. 

Scaling factors - Geometric similarity is ensured 
by scaling all dimensions in the ratio of the blade 
pitch at the tip, between the fan and the present inducer 
as follows: 
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Suction surface 



Figure 18 — Definition of Inducer Geometric Parameters 
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s = pitch =■ 


2tt r 


It 


where : 


= leading edge tip radius 

Z. = blade number 
1 


„ 2tt 78.6206/2 , „ , 

S_ = — — - = 10.2914 cm 

fan 24 


(Reference 2, Table 7, Page 68) 


S = 2tt X 4.3716 _ , 

^inducer 19 

scaling factor = !iSii2e£ = 

®fan 


4457 cm 
0.1405 


A number of fan cylindrical sections are given in 
Appendix 9 of Reference 2. These sections were scaled by 
the above factor to provide the basic blade sections to 
be used in the present design, i.e.; 


(X,Y) 


inducer 


^inducer 

^fan 


{X,Y) 


fan 


where ; 

X,Y = coordinates of blade given in Appendix 9 

of Reference 2. These blade sections were 
"laid out" to 20x scale. 

Fan stagger angles - In Appendix 9 of Reference 2, 

the static fan stagger angles (C , . • ) each cylin- 

s 1 ic 

drical section are given. The rotating stagger angles 
?rot f'the change is due to the untwist of the blades 
caused by rotating inertia forces) were obtained by 

substracting an amount rotating* 

static stagger angles. CB,^ ^ratic " rotating* 
deduced from Figure 20 (Reference 2) . The basic 
assumption implied is that each section rotates 
as a solid body about the leading edge. 
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Fan A . values - A . is the "throat" dimension 
min mm 

for a cylindrical section. A . was measured from a 

min 

drawing wherein two adjacent fan blade profiles are 

drawn at the appropriate fan pitch (from Appendix 9 of 

Reference 2 and stagger angle as defined above) . 

A , then corresponds to the "throat" dimension 
min 

and was measured from the drawing. This was easily and 
accurately done by using the 20x section drawings for the 
present design. If is calculated using the same 

scale geometry, then the A . to A , ratio will be the 

same as for the fan. An important feature in the above 
approach is that the fan stagger angles were retained. 

Fan blade setting angles - The blade setting angles 
on cylindrical surfaces are needed for the fan rotor for 
calculating and the new stagger angles. Unfortu- 

nately, Reference 2 does not give the blade setting 
angles on cylindrical surfaces, but only on conical 
streamline surfaces used in the aerodynamic cal- 
culations. The cylindrical surface values were checked 
in an approximate way and were found to differ by at most 
1° from the conical surface values. Furthermore, the 
largest difference occurs at the least important location, 
that is at the hub where the approach relative Mach 
number is low. Consequently, the blade setting angle 

distribution (Figure 19) as given in Reference 2 was 
used to generate the A^^^^ values . 

Fan incidence angle distribution - The consistent 
fan incidence distribution with span is that corres- 
ponding with the design tip Mach number (1.2) and tip 
incidence angle (4“) . For the present design, the design 
tip Mach number of 1.2 corresponds to approximately 110% 
of design speed (N^^) for the progenitor fan. 

The incidence variation with % span for the present 
design tip conditions, taken from Reference 1, is shown 
in Figure 20. 
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% Span (from hub) 


Figure 19— BlodeMeanLine Angle Distribution of Rotor Leading 
Edge (on Stream Surfaces)- P8iW Blading 






incidence, i (degrees) 
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Fan A 


values - The approach A 


values were 


rel 

calculated from 


rel 


\el = ® cos 6 

where B is the fan section air angle calculated from 


S « 3 


b 


i . 

m 


is picked off Figure 19. 

The each section for the 20x scale geometry was 

calculated. The A . /A , ratios for the progenitor fan 

min rel 

then follow from the previously measured A . 's 

min 


Up to this point, all the required fan parameters 

have been generated. From this information and the 

design intent of maintaining constancy of these parameters, 

the quantities A . and ? ^ (which together with the 

min !iroti 

scaled profiles define the blading) were calculated as 
indicated below. 


Inducer A . values - The A . values for each 
min min 

cylindrical section are calculated from the constancy of 

A . /A , and the calculated value for A , . The latter 
min rel rel 

quantity is calculated from 

A T = s cos 3 
rel 

where S is the approach flow air angles deduced from the 
potential flow calculations and s is the blade pitch for 
the particular section. 

Stagger angles, “ The new section stagger 

angles were calculated in the following way. The blade 
mean line angle (at the inducer leading edge) was deduced 
using the potential flow air angles and the i^ incidence 

distribution from Figure 20. The change in blade nean 
line angle (between present design and the progenitor fan) 
was then added to the fan rotating stagger angle to pro- ■ 
duce the inducer rotating stagger angle. Implicit in 
this is that the blade sections are re-staggered about 
the leading edge. The relevant formula is: 
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^rot ^rot fan 


where : 

A 

with 6 k = S + 1 


fan 


m fan 


The airfoil shape derived by these design methods 
was used to define the inducer blade at 100% design speed. 
A radial element inducer design was used and this will 
prohibit deformation of the inducer blades between static 
and rotating conditions. Therefore the values can 

be used to specify the manufacturing coordinates of the 
blades . 


Blade layout - The blade profiles (for each section) 
for the present design were generated from the above 
information in the following way. 

The scaled blade profiles were set at the calculated 
?rot two suction surface profiles spaced 

one blade pitch s apart. 

The locus of A , values was then laid off from the 
min 

suction surface toward the leading edge of the "trailing" 
blade contour (see the sketch below) . 


Blade 



52 



The leading edge of the blade was constructed by 
joining the existing suction surface with a pressure 
surface contour at the leading edge satisfying the 

criteria. This specifies the leading edge normal thick- 
ness. 


Furthermore, the nominal blade mean line pressure 
and suction surface locations are known for each section 
and hence may be laid out. This information is available 
because the blade setting angle and blade thickness dis- 
tributions are prescribed together with the requirement 
that radial elements exist in the inducer. 

The final blade suction surface profile was gener- 
ated by blending from this prescribed suction surface 
into the original fan blade profile. The location of 
this "blend" point (tangent point) was carefully con- 
trolled near the tip regions of the blade where tran- 
sonic performance will be obtained. The location 

is approximately where the bow shock from the leading 
edge will strike the suction surface in transonic super- 
sonic flow. t^rcl approach ^1^ where A shocks vi?ill 
terminate on the suction surface, from imbedded super- 
sonic regions, in transonic subsonic flow approach 

< 1) . It was therefore important to keep the blend 
point well downstream of the location. 


Suction 

Surface 



Normal 

Thickness 


^Mean line 
P & W Blade Profile 
X 

New Pressure Surface Leading Edge 
Profile 
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Finally, the remaining pressure-surface contour was 
generated from the leading edge portion defined above to 
meet the specified blade pressure surface contour. 
Usually, this presents no great problem as a smooth con- 
tour can be so joined. 

For blade sections close to the hub, it is more 
difficult to match the auction surface contours of a 
radial element inducer blade to the original fan profiles 
except close to the leading edge. This is because the 
fan blade is not a radial element blade and has signifi- 
cant turning toward (or beyond) axial in the hub sections 
Nevertheless, it was possible to maintain the 

criteria for all section profiles. Furthermore, it is 
probably not critical to hold the exact blade profiles 
close to the hub since these are low Mach number blade 
sections . 

The final blade shape was obtained by laying out 
radial plane sections (sections normal to the rotor axis) 
of the blade and smoothly joining the pressure/suction 
surfaces through the original fan sections laid out as 
above. Some iteration is required to smooth out small 
irregularities that appeared in the radial plane profiles 
This was accomplished by slightly "shifting" some of the 
cylindrical section profiles. 
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AERODYNAMIC DESIGN OF IMPELLER 


Approach 

The aim of the impeller design was to generate the geo- 
metric details of the impeller blades and the hifs and cover 
meridional contours within the geometric constraints deter- 
mined in the Preliminary Design Phase for the 3 = - 30 * case. 

The inducer design actually forms part of the impeller design 
and is faired into the impeller form in the approach used here. 

The basic criterion used for the impeller design is to 
ensure that certain aerodynamic features are achieved that will 
yield the tip conditions set in the Preliminary Design Phase 
(PDP) . At the same time the blade loading distribution must be 
controlled. These aerodynamic features are influenced in some 
complex manner by the required geometrical quantities, i.e. the 
hub and cover meridional contours , the thickness distribution 
of the blades and the blade setting angle distribution. It is 
impossible to study all geometric combinations of these vari- 
ables. Therefore, certain geometric features (e.g. thickness 
distribution) are prescribed and maintained, while other 
features are altered to meet the aerodynamic criteria. 
Inevitably, a measure of arbitrariness enters the design when 
these fixed geometric features are prescribed. 

The available geometric choices may be limited by questions 
of impeller stresses, blade and impeller deflections, impeller 
inertia, blade vibrational modes and instrumentation constraints 
{such as have been impressed upon this design by the use of LDV 
instrumentation) . The aerodynamic criteria should be achieved 
within these constraints. This may lead to some compromises 
being made in the aerodynamic criteria or even in the basic 
specifications . 

For the present design, the use of LDV instrumentation 
played a major role in the design as discussed fully in the 
Introduction. The impact of the LDV constraints on selection 
of a minimum impeller tip depth and the no-splitter require- 
ment had already been considered in the Preliminary Design 
phase, as had the requirement for a constant inducer cover 
radius. As it transpired, these restrictions did not appear 
to prevent the attainment of the aerodynamic criteria. 
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The only methods presently available for the calculation 
of flow in the impeller are those assuming potential flow. 

The Creare flow model dictates that the flow in most of the 
impeller consists of an isentropic jet which, however, may 
separate. The separated jet will then have a constant Mach 
number on its separated face and a wake region (most likely 
of increas ' ,\g area) growing towards the tip. The potential 
flow solution must reflect this flow model. 

The remainder of this chapter contains a discussion of the 
potential flow program and the aerodynamic criteria to be used. 
The results of the calculations are then presented, together 
with the definitions of all the relevant geometry for the 
impeller. 


VANCO Potential Flow Program 

The potential flow program used is a version of the calcu- 
lation procedure of Vanco (Reference 8) , who extended the basic 
program of Katsanis (Reference 9) to deal with centrifugal com- 
pressor flows. Briefly, the program solves the velocity-gradi- 
ent equation with the assiomption of a hub-to-cover mean stream 
surface along arbitrary quasi-orthogonals (i.e. , lines which 
intersect each streamline once only) in the meridional plane. 

The velocity-gradient equation is Euler's equation transformed 
to a coordinate system fixed to the mean stream surface. The 
latter is the two-dimensional stream surface running up the 
middle of a blade channel and extending from hub to cover. For 
the VAMCO program, the quasi-orthogonals must be lines of con- 
stant 0 extending from hub to cover (they lie in the meridional 
plane by definition of 6 = constant) . The use of quasi-ortho- 
gonals has a clear advantage over methods which employ normals 
to streamlines in that they remain fixed straight lines regard- 
less of any streamline form change. Having found the mean- 
stream-surface relative velocities, velocities on the pressure 
and suction surface are found using the rapid approximate 
method of Stanitz and Prian (Reference 10) . 

The basic inputs to the program, other than the obvious ones 
of mass flow, rotational speed, number of blades, etc., are the 
end point coordinates of the quasi-orthogonals with the 
associated wrap angles (6) and the blade-normal thickness at 
hub and cover- The program linearly interpolates along the 
quasi-orthogonals to find the normal thickness at intermediate 
points. In fact, any normal thickness distribution may be input; 
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the linear distribution was chosen as being convenient. 
Since the quasi-orthogonals are fixed straight lines, the 
mean blade shape is completely defined by specifying wrap 
6 as a function of meridional distance m for the hub and 
cover streamlines. Furthermore, the pressure and suction 
surfaces are defined since the normal thickness is 
specified. The actual blade geometry for the present 
design will be defined in the above manner, i.e., in the 
"VANCO implied" way. 

Further important inputs to the program are the 
exit flow angle (set equal to the jet angle) and the 
possibility of including losses by means of a distribution 
of loss of relative total pressure with meridional dis- 
tance . 


Modifications to Reflect Creare Flow Model 

The VAMCO program cannot be used with the usual input 
if it is to reflect the Creare flow model of an isentropic 
jet flow and a growing wake region. However, the basic 
physics of this model can be represented approximately by 
an appropriate change in the values of the normal input 
quantities as outlined below. 

For the present design, a fairly simple and tradi- 
tional approach was taken. No experimental evidence is 
available that indicates how the wake region grows within 
the blade passages of the centrifugal impeller. There- 
fore, recourse must be made to a cruder representation. 

For the present design, this entailed iterating boundary 
layer calculations with the potential flow calculations 
and then treating the boundary layer displacement thick- 
ness as wake flow by taking it as additional blade 
blockage . 


Aerodynamic Criteria 

The aim of the aerodynamic criteria is to ensure, 
insofar as possible, that the Preliminary Design Phase 
(PDF) tip conditions are matched. These rules are; 

1) Achieve the PDP target impeller tip diffusion 
ratio (DR 2 ) . 

2) Limit the diffusion to no more than 50% in 
the inducer (i.e, , near-axial portion of the 
impeller) . 
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3) Limit the diffusion to no more than 15% futher 
diffusion in the "radial" portion of the impeller. 

4) Ensure that, if separation does occur, it does so 
two-dimensionally . 

Recall that the traditional Creare flow model with the 

assumption of a constant relative Mach number along the jet face, 

starting at the jet/wake separation location, means DR_ < DR 

2 sep 

Clearly maximizing DR implies that DR^ is maximized. The 

former was the intent of the Preliminary Design Phase and so 
presumably Rule 1 is already accomplished, assuming that the 
remaining three rules are satisfied. 


The other rules seek either to ensure that separation occurs 
along a two-dimensional line between hub and cover suction 
surface or to eliminate separation altogether. Rule 2 suggests 
that no more than 50% diffusion be allowed in the inducer portion 
of the impeller or else the turbulent boundary layers will 
separate there. This separation cannot be accomplished two- 
dimensionally because the hub is far too low. The separa- 

tion line should be a line of constant M for steady internal 


Some recent work (Ref. 13) has shown that fully developed 
turbulent boundary layers in rotating channels may become 
"stabilized" when a critical rotation number is exceeded. 

"Stabilize" means the boundary layers take on a laminar 
character and it is postulated that no more than 15% diffusion 
may be applied before separation occurs. Actually separation 
is allowable in the radial portion because this can usually 
be done two-dimensionally . However, such a separation will be 
acceptable only if the separation Mach number (on suction 
surface) matches the target tip value from the Preliminary 
Design Phase. 

The most critical stream tube affecting the flow in the 
impeller is the one along the blade suction surface at the cover. 
When the above rules are applied to this stream tube, the loading 
distribution suggested is that of a rapid diffusion in the 
inducer or axial portion, follov;ed by a gradual diffusion in the 
radial portion. At the same time, the rules are applied to the 
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hub suction surface, so that no more than 50% diffusion 
should occur in the axial portion of the hub stream tube. 
Hub separation, should it occur, will tend to raise the 
cover suction surface relative velocity because of 
blockage of the separation bubble, making it more diffi- 
cult for the flow on this latter surface to reach the 
target M within the 15% radial diffusion constraint. 


The potential flow calculations for the pressure 
surfaces will show higher diffusion rates than for the 
suction surfaces, indicating the possibility of separa- 
tion on these surfaces. However, this will probably not 
occur since the pressure surface boundary layers are 
vigorously stirred in the radial portion of the wheel by 
the action of Coriolis forces, which remove low energy 
boundary layer fluid created on the blade driving surface. 
In the axial part of the impeller the driving mechanism 
is the large blade-to-blade loading. 


Another important constraint is that no "potential 
flow reversal" should occur. This is most likely to 
occur upon the hub pressure surface. Reverse flow at the 
hub has the effect (because of its pseudo blockage effect) 
of accelerating the flow on the suction surface. For 
these reasons, rapid deceleration at the hub, if any, 
should be delayed as far as possible toward the impeller 
exit. 


Creare's ideal loading distribution is indicated in 
Figure 21. For the present design, an attempt was made to 
match this loading distribution and to meet the target 
^rel 2j same time. If possible, the target 2j 

is to be met without separation occurring. If separation 
does occur, it should occur two-dimensionally. As may be 
seen, this can only be accomplished near the impeller exit 
if the ideal loading distributions are to be followed. 

These criteria were used in the present design. 
However, they were impacted by the mechanical constraints 
of hub and blade stress limitations and by the desire to 
have a constant radius at the inducer shroud where the 
LDV window is located. A low number of blades at 
impeller exit (i.e. no splitters) was desired for the 
LDV studies. 
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Results 


The above general approach was applied to the design of the 
present impeller. The geometric constraints, already fixed in 
the Preliminary Design Phase, are the inducer leading-edge hub 
and cover radii and , impeller tip radius, depth and 

blade angle (r_, b_ and B, and a constant inducer cover radius 
^ 2 . 02 

deep into the impeller (LDV constraint) . 

Within the above constraints , the following geometric 
features must be chosen: 

i) the hub and cover meridional contours, 

ii) the hub and co'^rer blade setting angle (3 and 3j^^) 

distributions, and 

iii) the hxib and cover normal blade thickness distributions . 

Initiation - The above items were set initially in the 
following way. The cover in the inducer portion was set at 
constant radius (i.e., order to eliminate the need for 

a curved inducer LDV window. The hub contour was made a 7° 
ramp in order to satisfy (approximately) continuity at the end 
of the inducer. The remaining part of the contours (both hub 
and cover to the tip) were generated using a simple continuity 
check. 

The cover 3 distribution was set by drawing a smooth curve 

from the known tip value of -30° (3.,, ) to blend into the 

2b 

previously prescribed inducer (3j^) distribution. In the inducer 

section, the hub 3. distribution follows from the radial element 

b 

assumption, while the rest is generated by drawing a smooth 
curve between the inducer portion and the tip value of -30°. 

The blade normal thickness along the cover was set in the 
inducer portion by laying out the tip inducer section and then 
using a constant blade thickness distribution for the remainder 
of the blade. The blade normal thickness distribution along the 
hub was initially chosen "arbitrarily" with a peak value (chosen 
by intuition) at the position likely to have the hiahest stress 
value, followed by a fall-off in thickness. 
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Iteration. - The basic procedure was to vary the shape of 
the meridional contours and the distributions in attempting 

to match the ideal loading distribution (see Figure 21) • The 
normal thickness distributions remained essentially fixed. The 
hub normal thickness was changed slightly to alleviate stress 
problems as these stresses were calculated. The final blade 
normal thickness distributions at hub and cover are shown in 
Figure 22. 

The blade turning distribution (3j^) has a strong effect on 

the loading distribution in the axial part of the impeller (i.e., 
inducer region and slightly beyond) , but less so in the radial 
portion where centrifugal and Coriolis effects largely determine 
the loading. Consequently, the distribution in the latter 

half of the inducer portion (on the cover side from Z = 0.89 cm 
to about approximately Z = 3.3 cm) is varied in conjunction with 
the meridional contours to produce the required cover suction 
surface loading distribution. Recall that the distribution 

in the first portion of the inducer is determined by the fan 
scaling (see Inducer Section) . 

The rest of the 6j^ distribution is developed by laying out 

a smooth curve from the tip value (-30°) to blend into the 
desired distribution (indicated above) just beyond the end 

of the inducer. This curve is not drawn completely arbitrarily 
in that it is necessary to keep the B ’ s and hence the blade tip 

wrap 0^ as low as possible in order to minimize blade stress. 

The hub distribution is largely determined by the cover 3j^’s 

because : 

1) radial elements are required in the inducer portion, 

2) zero rake angle (0 plane) at the tip is prescribed. 

This means that the hub blade wrap angle (0^ - Qj^h 

must be equal to the cover blade wrap angle (B^ - 

The final blade 3 distributions at cover and hub are shown 
in Figure 23. 
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Figure 23 — Blode Angie Distribution — Hub and Cover 



The meridional contours (particularly the cover contour) 
have a large effect on the mean and blade surface velocities. 

The principal effect is in the regions beyond the inducer 
because of the significant curvatures there. (This will be 
especially true in a low inertia impeller.) Initially, a 
reacceleration on the cover suction surface was encountered 
just beyond the inducer end. This was eliminated by increasing 
the impeller length, thereby decreasing the curvatures. Recall 
that the cover contour was constrained by the requirement of 
constant radius in the inducer section. The hub meridional 
contour is largely determined once the cover one is set because : 

i) the tip radius and depth (r_, b ) are known, and 

ii) the hub ramp angle and length are given. 

Boundary layer calculations - Boundary layer calculations 
were carried out for the cover and blade suction surfaces using 
a boundary layer prediction technique based upon the methods 
developed by Englert (Reference 11) . The input for these calcu- 
culations came from the potential flow results. The assumption 
was made that the boundary layer growth throughout the impeller 
is typical of that calculated for the blade suction surfaces. 

The action of Coriolis forces is to sweep excess boundary layer 
fluid created on the blade driving (pressure) surfaces toward 
the. suction surface near the exit of the impeller. 

Final results - As well as iterating the geometrical 
quantities indicated above in attempting to reach the ideal 
loading distribution, the potential flow calculations were 
iterated with the boundary layer calculations. This was accom- 
plished in the traditional manner of calculating the potential 
flow without boundary layers as a first step, then calculating 
the boundary layers as indicated above and re-doing the potential 
flow calculations by including the boundary displacement thickness 
as extra blade blockage. Only one loop of this iteration was 
performed. 

The final results are displayed in Figures 24 to 31 . 

Figures 24 through 29 show the relative Mach number and relative 
velocity distributions as calculated for the impeller flow. 
Distributions are shown as a function of impeller meridional 
distance for the cover streamline (streamline 13) , for the 
approximately mid-span streamline (streamline 7) , and the hub 
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Fiqure28- Relative Velocity Distribution — Mean Span Streamline 





streamline {streamline 1) . For each streamline, the data 
is displayed for the suction (trailing) blade surface, the 
pressure (leading) blade surface, and for the flow passage 
mean streamline. 

The calculations are shown "dashed" immediately at the 
entrance to the impeller, since the calculated results are 
not expected to be correct in the leading edge portion of 
the inducer blades. 

A rapid diffusion occurs in the axial inducer, 
followed by the more gradual diffusion in the radial portion 
of the impeller. In these figures, the cover and mid-span 
plots show an initial, large "loading" in the inducer and 
then a reduced loading until the loading is again increased 
in the radial portion of the impeller. Actually, the 
loading is expected to be more uniform, e.g., as shown 
in Figure 24; the spline fit routines used in the calcu- 
lations are incapable of accurately following the details 
of the blade g distributions in this region. Thus, the 
ideal loading distribution (Figure 21) was matched reason- 
ably well for the hub streamline, but the cover streamline 
exhibits too much loading in the radial part of the impeller. 
The use of splitters would have prevented the latter. 

Separation is expected to occur on the cover suction 
surface near the exit of the impeller at the target separa- 
tion relative Mach number of 0.792. The flow is assumed 
to establish a separation line across the blade span at 
this Mach number and to flow at constant jet Mach number 
to the impeller exit. 

There is a significant divergence between the actual 
and ideal loading near the impeller tip. However, the 
actual loading, velocity, and Mach number distribitions 
at and near the impeller exit are probably not correct 
because (1) the flow is expected to separate within the 
impeller as discussed above, and (2) the blade loading 
calculation routine calculates loading at the impeller 
exit where none will exist and the routine is "sensitive" 
to rapid changes in flow blockage (such as exists near 
the impeller exit when the boundary layer displacement 
thickness is included) . 

Figures 30 and 31 show the cover and hub blade 
boundary layer displacement thickness and shape factor 
H (H 6*/0) distributions. These distributions have been 
used for both pressure and suction surfaces in setting 
the effective boundary layer blockage distribution used in 
the potential flow calculations. A linear distribution 
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DisplQceme 











between hub and cover values was used across the blade span. 
Figure 30 shows a rapid increase and decrease in boundary 
layer characteristics in the middle of the inducer on the 
cover streamline. The shape factor far exceeds traditional 
boundary layer separation values (of order 2.7). However, 
as indicated above, the rate of change of blade angle (and 
hence diffusion rate) will be less than calculated and, 
hence it is anticipated that separation will not occur in 
the inducer. 


Final Blade and Passage Definition 

"Nominal" hub and cover meridional contours are defined 
corresponding to the above calculations. Additional calcu- 
lations must be made before the final coordinates for the cover, 
the blade tip, and the hub contours are defined. The coordi- 
nates for the "nominal" hub and cover meridional contours are 
given in Table VI together with the blade definition (8, t^^) at 
the same points. 

Turbulent boundary layer calculations were carried out 
for the surfaces of the compressor inlet (both hub and cover) 
using the Englert (Ref. 11) method (see below for the inlet 
contour definition) . These calculations indicated approxi- 
mately the same displacement thickness at the inducer leading 
edge for both the outer and inner inlet contours. The 
boundary layer blockage effect must be considered in defining 
the cover and hub meridional contours. 

Cover contour - The final cover contour was constructed 
very simply by using the nominal cover contour as a base. A 
smooth contour was constructed from a point displaced radially 
0.025 cm at the inducer leading edge to bland into the nominal 
cover contour just beyond the end of the inducer. This curve 
defines the final cover contour. 

The blade tip contour is defined simply by specifying a 
clearance distribution with respect to the final cover contour 
as defined above, A 0.025 cm running clearance normal to the 
above meridional contour was used from impeller inlet to exit. 
This represents a clearance of approximately 8% of the impeller 
tip exit depth. 

Hub contour - The actual hub contour was constructed by 
reducing the leading edge radius by 0.023 cm (to allow for 
boundary layer blockage) and putting in a ramp angle of slightly 
more than 7® to blend into the nominal hub contour just beyond 
the inducer end. Figure 32 shows the impeller meridional view. 

Blade geometry - For the inducer portion of the impeller, 
the blading was developed by scaling the P & W fan blading as 
stated in the Inducer Section. 
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TABLE VI 

: - NOMINAL IMPELLER GEOMETRIC ' 

COORDINATES^ 



Quasi- 

Blade 

Axial 

Radial 

Angular 

Normal 

Orthogonal 

Angle , 

3 

Distance, X 

Distance , r 

Location , 0 

Thickness, t 
, , n 

Number 

(degrees) 

(cms) 


(cins j 


(degrees) 

(cms) 


cover 

hub 

cover 

hub 

cover 

hub 


cover 

hub 

4 (inlet) 

-56.34 

-36.70 

0. 

0. 

4.3716 

2.1858 

0. 

0.0152 

0.0305 

5 

-54,24 

-35.33 

0.5080 

0.5080 

4.3716 

2.2482 

-9.2537 

0.0610 

0.1829 

6 

-43.33 

-26.33 

1.0160 

1. 0160 

4.3716 

2.3105 

-16.9675 

0.0610 

0.2489 

7 

-37.18 

-22.36 

1.6510 

1.6510 

4.3716 

2.3885 

-24.0015 

0.0610 

0.3048 

8 

-30.25 

-19.56 

2.1755 

2.1538 

4.3859 

2.4603 

-28.6043 

0.0610 

0.3401 

9 

-24.78 

-17.33 

2.6708 

2.5853 

4.4294 

2.5440 

-31.9627 

0.0610 

0.3602 

10 

-19.94 

-14.94 

3.2841 

3.0812 

4.5345 

2.6691 

-35.2184 

0.0610 

0.3759 

11 

-17.21 

-12.97 

3.8417 

3.5220 

4.6954 

2.8124 

-37.6289 

0.0610 

0.3833 

12 

-16.03 

-11.21 

4.-3413 

3.9443 

4.9142 

2.9855 

-39.5633 

0.0610 

0.3912 

13 

-16.11 

- 9.06 

4.8808 

4.5064 

5.2727 

3.2892 

-41.6539 

0.0610 

0.4013 

14 

-17.14 

- 7.24 

5.2419 

5.0461 

5.6300 

3.6971 

-43.2458 

0.0610 

0.4094 

15 

-18.40 

- 5.99 

5.4701 

5.5383 

5.9407 

4.2394 

-44.4685 

0.0610 

0.4163 

16 

-20.06 

- 6.76 

5.6607 

6.0262 

6.2889 

5.1405 

-45.7654 

0.0610 

0.4191 

17 

-22.13 

-10.10 

5.8128 

6.2802 

6.6756 

6.1284 

-47.1816 

0.0510 

0.4155 

18 

-25.55 

-18.85 

5.9541 

6-3500 

7.2527 

7.1275 

-49.3400 

0.0610 

0.3353 

19 

-30.00 

-30.00 

6.0262 

6.3500 

7.9756 

7.9756 

-52.2312 

0.0610 

0.2032 


^Before correction on hub and shroud for inlet boundary layer growth. 






Figure 32- Meridional View of Impeller 


The blading for the remainder of the impeller was develt- 
oped in the "VANCO implied" way. Thus, radial plane cuts 
were passed through innumerable quasi-orthogonals (constant 
0 lines) , which define the mean line of the blades. The 
pressure and suction surfaces were then developed assuming 
a linear distribution in normal thickness along the quasi- 
orthogonals with a linear blade angle distribution along the 
radial plane cuts. Some smoothing of the blade contours was 
done by iterating between tangential plane cuts and the 
radial plane cuts. 


Inlet Design 

As stated above, boundary layer calculations showed 
approximately a 0.025 cm boundary displacement thickness on the 
outer and inner inlet contours at the inducer leading edge. 

This represents an inlet area blockage of approximately 2.25% 
based on flow area at the impeller inlet. This blockage was 
taken into account in defining the inlet contours shown in 
Figure 33. The final inlet configuration uses long radius 
ellipses on both the inlet "bellmouth" outer contour and the 
inner centerbody contour. The centerbody is supported on an 
inlet sting of 0.935 cm diameter extending from upstream in 
the inlet plenum chamber. Both the bellmouth and centerbody 
configurations have been used in previous compressor test 
programs. Adequate aerodynamic and mechanical performance 
has been obtained with these designs. 
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DIFFUSER FINAL DESIGN 


Objectives 


The objectives of the final diffuser design exer- 
cise were: 

1) to establish the diffuser geometry in detail, 

2) to identify any deviations from the prelim- 
inary design phase performance predictions 
which are occasioned by the final diffuser 
design geometry, and 

3) to examine and evaluate an alternate diffuser 
design which could lead to possible improved 
diffuser performance in terms of lower losses 
and increased range. 

Input from Preliminary Design Analysis 

All of the geometric specifications for the dif- 
fuser, as established by the preliminary design phase 
output were used in the design of the primary diffuser 
configuration. The second diffuser configuration, called 
here the alternate design configuration, used an in- 
creased diffuser depth. This increased depth, with the 
constraint of maintaining constant throat Mach number and 
hence total throat area at the design point, resulted in 
an altered channel diffuser throat aspect ratio. Hence 
a re-evaluation of the channel diffuser geometry was also 
undertaken. 

The primary diffuser design had as basic input from 
the Preliminary Design Phase all the geometry and aero- 
dynamic quantities for the B_, = -30° case (see Table V) . 

2b 

The diffuser vane-leading edge radius was not estab- 
lished definitively by the preliminary design analysis. 

A radius ratio - 1.05 was considered to be approxi- 
mately that desired based upon Creare experience. The 
choice cannot be defended on a very sound deterministic 
fluid dynamic basis. Small values of are known to: 

1) increase noise level, 

2) increase vibratory stress of impeller-blade 
trailing edges. 
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3) increase vibratory stress in the diffuser- 
vane leading edge region, and 

4) increase the strength of the aerodynamic 
generated pressure fields in the impeller 
exit/diffuser entry region. 

The limit on increasing the ratio is basically 

set by wall friction. While the design correlations for 
stagnation pressure loss between impeller tip and 

channel throat, and impeller backflow loss do not now 
recognize (nor any other diffuser geometrical values) , 

it is expected that tends to rise above the blockage 

correlation as increases. In general, the assumption 

is made that it is undesirable to increase the extent of 
the sidev/alls of the vaneless/semi-vaneless region; this 
implies keeping R^ as small as possible. 

As will be seen in the discussions that follov/, the 
method used for defining the diffuser vane suction sur- 
face has effectively set the value of R_ at 1.0764 
(r^ = 8 . 585 cm) . ^ 


Approach 


Primary diffuser design - The primary diffuser ge- 
ometry was generated by applying the fluid dynamic model 
and methods described below to the diffuser entry region 
at 100% Np and at design mass flow. It is to be recog- 
nized that the diffuser geometry generated by this 
approach, while presumably satisfying the aerodynamic 
performance at the design point conditions, may not nece- 
ssarily represent the optimum geometry for other oper- 
ating flows and impeller speeds. For example, incidence 
of the flow onto the diffuser vane lead edge, a change of 
channel diffuser throat Mach number with flow, etc. will 
change the level of overall diffuser recovery and hence 
stage performance. 

The basic design approach assumes a core potential 
flow in the vaneless and semi -vaneless regions with tur- 
bulent boundary layer development on the cover and hub 
diffuser walls and on the vane suction surfaces. An 
iteration approach was used between successive solutions 
of the boundary ^ layers and the core potential flow to 
arrive at the flow solution for the vaneless and semi- 
vaneless regions of the diffuser. 
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Once the boundary layer/core flow solution was 
obtained, the vane tip radius and vane suction surface 
were prescribed. These were set to provide for a smooth 
deceleration of the flow along the vane suction surface 
and a flow Mach number distribution approaching the 
channel diffuser throat close to that prescribed for the 
one-dimensional throat Mach number from the preliminary 
design analysis. 

The potential core flow was calculated as a com- 
pressible, axisymmetric vortex flow. The flow was 
assumed completely mixed out to the 2 * state at the 
impeller exit. The compressible vortex solution satis- 
fies the mass flow and angular momentum relations at any 
radius ; 

mass flow = 2irpC (b-26*)r 

in 

angular momentum = 2 iTr^pC^Cg (b- 2 fi*) 

where : 

Cjjj and C 0 are the local core flow radial and 

tangential velocities respectively, 
b is the local diffuser depth (between cover 
and hub walls) , and 

6 * is the boundary layer displacement thickness. 

The Englert (Reference 11) compressible, turbulent 
boundary layer calculation technique was used to calcu- 
late the boundary layers on the hub and cover walls. The 
initial conditions for the boundary layer calculations 
were as follows. A turbulent boundary layer was pre- 
scribed to "start" on the hub wall since the boundary 
layer fluid leaving the impeller is "thrown" onto the 
hub wall at essentially wheel speed. The cover wall 
boundary layer was assumed to have a thickness equal to 
the prescribed running clearance between the impeller and 
cover at the impeller tip; it was assumed that the bound- 
ary layer cannot be any thicker than this since the 
impeller blades effectively "chop off" the boundary layer. 

An empirical correlation (Ref. 1) of throat blockage 

versus non-dimensional pressure rise from the impeller 
exit-to-dif fuser throat, used to establish 

the magnitude of boundary layer blockage at the throat. 

The boundary layer calculations, including the boundary 
layer blockage calculated on the vane suction surface, 
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usually do not agree exactly with the empirical vs 
Cp 2*_4 correlation. Therefore, the boundary layer displace- 
ment thickness distribution between the impeller exit and 
channel throat was varied, proportional to streamline dis- 
tance from the impeller exit, by an amount sufficient to 
produce the throat blockage given by the empirical correla- 
tion. Thus the throat area, aspect ratio and mean Mach 
number were all maintained as given by the preliminary 
design calculations. 

The vane suction surface was displaced away from the 
core flow "suction surface streamline" by an amount equal 
to the calculated displacement thickness of the suction sur- 
face boundary layer so as to maintain the core flow stream- 
lines as calculated. Finally, the location of the vane tip 
radius was prescribed by the necessity of the tip radius 
and vane suction surface contour providing the required 
channel diffuser throat width, W^. The vane tip radius was 

set as small as practical from fabrication and assembly 
considerations. The vane suction surface (plus boundary 
layer displacement thickness) is aligned to follow the 
vortex flow streamline. The vane pressure surface was laid 
out tangent to the vortex flow streamline impinging on the 
vane lead edge. Thus all the divergence is taken in the 
suction surface. The mass averaged Mach number at the 
diffuser exit was calculated to be 0.155 with a total 
pressure of 8.24 x lO^N/m^. 

Figure 34 displays the type of vaned diffuser geometry 
specified by this design procedure and summarizes the condi- 
tions established by this design approach. 

Alternate diffuser design - An alternate diffuser con- 
fig uraOTon~sepixate~~froiir~the~priraary diffuser design, was 
designed and fabricated. This diffuser is called the step 
diffuser because a sudden expansion (increase in axial 
diffuser depth) is provided immediately at the impeller exit 
radius. The concept behind this diffuser design is the 
provision of a "smoke shelf" on the front diffuser cover at 
the impeller exit radius. This "smoke shelf" is intended to 
provide the effect established in the chimney of a fireplace, 
namely, to establish a fixed separation point for the flow 
and prevent the backflow of fluid. In this case, the "smoke 
shelf" or step is provided on the diffuser front cover to 
prevent radially momentum deficient fluid on the cover from 
being driven back into the impeller by the imposed diffuser 
radial pressure gradient. 
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D i ffuser porometers 



f 2 = 79.76 
Zq= 17 
rj = 85.85 
tC 3 = 0.254 

f s o*> 

R = 70.04® 

^b3 

W 4 * I 1.48 
64 = 3.24 
29= 10 ® 

L/W 4 = 17 
rg = 254 
0 21. 18® 

Circumf erentiol 
plane divergence 

All dimensions 
in mm. 


Figure 34— Voned D iff user Geome tr y — Pr i mory Design 






The design provided for a step height equal to 10% of 
the nominal impeller exit depth. Downstream of this step 
change in depth, the diffuser passage was held at constant 
depth. The remaining design procedures, namely the 
specification of the diffuser vane suction surface geometry, 
lead edge radius and channel diffuser geometry, followed 
the same design methodology as described above for the 
primary diffuser. The additional input required in the 
step diffuser design was the description of the boundary 
layer displacement thickness growth on the front cover in 
the region of the step separation. Correlations of separa- 
tion length behind backward facing steps, under a range of 
free stream Mach numbers up to sonic, show that the separa- 
tion length is on the order of 5 to 10 step heights (Ref. 

14 to 15). For the step height of this design (3.23 x 10“"* 
meters) , the length of the separation zone is so small, order 
of 3.23 X 10“^ meters, that the influence of the free shear 
layer on the growth of the displacement thickness over the 
separation zone has been neglected; the boundary layer has 
been assumed to grow in displacement thickness the same as 
would occur if the separation zone were a solid surface. 

The final diffuser vane geometry was changed from the primary 
diffuser design because the free vortex flow occurs with a 
wider diffuser depth. 

Figure 35 displays the final vaned diffuser geometry 
specified by this design procedure and summarizes the condi- 
tions established by this design approach. 
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Diffuser porometers 



f2= 79.76 
2 q= I 7 

85.48 
tC3= 0.Z54 
0“ 

= 7i.87« 

b3 

%= 10.43 
b4= 3.56 
20 = 10 "* 

L/%= 18.6 
rs = 254 
e^= 21. IS” 

Circumferential 
plane divergence 

All dimensions in mm. 


Figure 35-Voned Diffuser Geometry —Alternate Design 


COMPRESSOR INSTRUMENTATION 


The compressor unit has been instrumented (emphasis 
on aerodynamic instrumentation) in order to provide (i) 
an adequate description of the flow states within the 
compressor unit and (ii) monitor shaft vibrational modes 
during running. This instrumentation includes instanta- 
neous and time-averaged static pressure taps, time- 
averaged total pressure probes, temperature thermocouples, 
inlet flow survey instrumentation and dynamic clearance 
probes. These items are fully described below (with 
figures and tables) according to type of instriment and 
location. 

Also included in the final instrumentation package 
to be used on this compressor is a Laser Doppler Veloci- 
meter (LDV) system (designed and fabricated under a sepa- 
rate program) . This LDV system will be used to measure 
velocity components in the impeller inducer region and in 
the diffuser vaneless and semi-vaneless spaces. Special 
LDV windows have been inserted in the cover and in the 
diffuser front plate to provide quality optical surfaces 
for the LDV optical studies. The critical LDV window 
configuration is the diffuser one, where the windov; must 
span several diffuser vane tips and mechanically support 
and aerodynamically seal across the diffuser vanes. 
Analyses and experimental bench tests have been made of 
these LDV window configurations to ensure that the prob- 
lems associated with the thermal mechanical stress, dif- 
ferential thermal expansion and resultant leakage between 
window quartz and diffuser cover metal, and mechanical 
clamping of diffuser vanes with the quartz window can be 
successfully dealt with in compressor test operation. 

The design of the LDV system has been carried out under 
a companion contract. Full details of the LDV configu- 
ration may be found in Reference 12. 

Type and Location of Aerodynamic Instrumentation 

The following measurements will be made on this 
compressor unit: 

- time-averaged static pressures 

- instantaneous (time varying) static pressures 

- time-averaged total pressures 



- temperature measurements 

- inlet flow direction survey 

Reference system - Top dead center on the machine 
is identified as 0 = 0°. Theta is measured positive 
counterclockwise when looking into the inlet. One dif- 
fuser vane tip is located at 0 = 0“- An odd number of 
vanes then implies that bottom dead center splits the 
pitch of two vane tips. Since pitch = 21. 18", 0 = 0° 
aligns with the first vane tip and 0 = 10.6“ is mid-pitch. 
The inducer blade lead edge is set at X - 0, and X is 
positive in the flow direction. 

Time-averaged static pressures 

a) Inlet Static Taps - See Figure 36. 

Three (3) taps, each at two (2) circumferen- 
tial locations were positioned on the inlet 
cover. The axial locations of these taps 
are 0.63, 1.27, and 2.54 centimeters ahead of 
of the inducer leading edge (X = 0) . The 
circumferential locations are at 6 = 0“ and 

e = 10.6“. 

b) Impeller Cover Taps - See Figure 37. 

Impeller cover taps are located at three (3) 
circumferential locations. These are at mid- 
pitch of channel 1 (0 = 10.6°), 5/6 pitch of 
channel 1 (0 = 17.66“), and 1/6 pitch of 
channel 2 (0 = 24.72“) . 

The meridional locations of these taps 
coincide with the shroud streamline constant 
0 locauions used in the VANCO program, plus 
several extra tap locations. Table VII pre- 
sents the locations of these taps on the 
cover inside surface in terms of X (axial 
distance from inducer leading edge) and 
r (radius) . 

c) Impeller Tip Static Taps - See Figure 38. 

Impeller tip static pressure taps are located 
on the front cover in an arc to span a diffuser 
entry region. Eight (8) taps span the entry 
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Figures?— ImpeHer Cover Stotic Tap Locations 










TABLE VII - IMPELLER COVER STATIC TAP LOCATIONS 


Tap Number 

Axial 

Distance 

(cms) 

, X Radial Distance, r 

(cms) 

7a 

7 

22 


0 

4.40 

8a 

8 

23 


0.51 

4.40 

9a 

9 

24 


1.02 

4.40 

10 a 

10 

25 


1.65 

4.40 

11a 

11 

26 


2.18 

4.41 

12a 

12 

27 


2.67 

4.45 

13a 

13 

28 


3.28 

4.56 

14a 

14 

29 


3.84 

4.72 

15a 

15 

30 


4.34 

4.94 

16a 

16 

31 


4.88 

5.30 

17a 

17 

32 


5.24 

5.65 

18a 

18 

33 


5.47 

5.97 

19 a 

19 

34 


5.66 

6.31 

20a 

20 

35 


5.81 

6.70 

21a 

21 

36 


5.95 

7.25 

Taps 

7a 

through 21a are 

located 

at 0 = 10.6°. 






See 

Taps 

7 

through 

21 are 

located 

at 6 = 17.66°. Figure 37 


Taps 22 through 36 are located at 9 = 24.72° 




Figure38~ Pressure Static Taps and Semi-Conductor Transducer 


Locations in Diffuser inlet 



region between vane tips (in channel 1) . 
Another set of eight (8) taps has been placed 
at r = 1.02 radius. This second set is 

located one vane island pitch over from the 
first set at r = 1.0 r2 (i.e., in channel 2). 

r, 9 positions are given in Table VIII. 


TABLE 

VIII - IMPELLER TIP 

STATIC TAP LOCATION 

Tap Number 

Radial Distance, 

r Angular Location, 6 

(cms) 

(degrees) 


37 

7.99 

0 

38 

7.99 

2.71 

39 

7.99 

5.57 

40 

7.99 

8,24 

41 

7.99 

10,91 

42 

7.99 

13.40 

43 

7.99 

15.97 

44 

7.99 

18.74 

45 

8.16 

22.70 

46 

8.15 

25.34 

47 

8.15 

28.05 

48 

8.15 

30.86 

49 

8.14 

33.55 

50 

8,15 

36.07 

51 

8.15 

38.58 

52 

8.15 

41.12 


d) Vaneless/Semi-Vaneless Array - See Figure 38. 
An array of static pressure taps is located 
in the vaneless and serai- vaneless region 
between the impeller tip and the diffuser 
throat on the front cover. The taps form 
three (3) arcs which approximate the suction 
surface, centerline and pressure surface 
streamlines for channel 1 . Table IX dis- 
plays the r, 9 coordinates of these taps. 
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TABLE IX - VANELESS/SEMI-VANELESS STATIC TAP LOCATIONS 


Tap Number 

Radial Distance, r 
(cms) 

Angular Location, 6 
(degrees) 

53 

8.51 

-1.00 

54 

8.64 

1.27 

55 

8.76 

3.30 

56 

8.87 

5.54 

57 

8.99 

7.76 

58 

9.12 

9.91 

59 

9.25 

12.16 

60 

9.36 

14.27 

61 

9.48 

16.45 

62 

9.63 

18.43 

63 

9.80 

20.23 

64 

9.98 

22.05 

65 

8.20 

2.20 

66 

8.32 

4.30 

67 

8.44 

6.66 

68 

8.56 

8.99 

69 

8.68 

10.99 

70 

9.25 

13.13 

71 

8.91 

15.13 

72 

9.02 

17.32 

73 

9.16 

19.59 

74 

9,28 

21.67 

75 

9.43 

23.62 

76 

8.21 

12.24 

77 

8.33 

14.38 

78 

8.45 

16.57 

79 

8.55 

18.67 

80 

8.69 

20.96 

81 

8.80 

23.04 

82 

8.95 

25,11 

83 

8.63 

-0.40 

102 

8.55 

42.50 


e) Channel Diffuser Taps 

An array of static pressure taps is located in 
one channel diffuser on the front cover. These 
taps are spaced approximately three (3) throat 
widths apart along the centerline of channel 1, 
Table X presents the r, 6 coordinates of the 
taps as well as those of additional throat taps. 
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TABLE X - 

■ CHANNEL DIFFUSER 

STATIC TAP LOCATIONS 

Tap Number 

Radial Distance 

, r Angular Location, 0 


(cms) 

(degrees) 

84 

10.81 

35.51 

85 

13.16 

46.85 

86 

15.80 

54.41 

87 

18.60 

59.72 

88 

21.73 

63.77 

89 

24.72 

66.69 

90 

9.18 

-1.98 

91 

9.17 

40.77 

92 

9.17 

61.97 

93 

9.16 

83.22 

94 

9.17 

104.60 

95 

9 . 15 

125.78 

96 

9.16 

147.00 

97 

9.17 

168.20 

98 

9.17 

230.90 

99 

9.17 

252.14 

100 

9.18 

266.76 

101 

9.20 

336.90 

Time Varying 

Static Taps 



All dynamic (instantaneous) pressures are measured 
with semi-conductor pressure transducers. These trans- 
ducers have a 0.216 cm diameter pressure-active diaphragm 
with an allowable minimum spacing between transducers 
of 0.939 cmr center to center. Because of a large zero 
pressure shift with temperature, DC pressure levels are 
not easily determined with these transducers. 

a) Inlet and Impeller Locations - See Figure 39. 

There are ten (10) semi-conductor pressure 
transducer locations spaced through the 
inducer and impeller. The scheme employed 
was to wrap the transducer locations to 
correspond with the impeller blade wrap so that 
the transducers see the blades beneath them 
simultaneously. The locations are staggered 
circumferentially as required for clearance 
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F- igure39~lnlet and Impeller Sem i- Conductor TronsducerLocotions 





and so as not to interfere with the time- 
averaged static taps. One probe is located 
0.635 cms ahead of the inducer leading edge. 
The coordinates of these taps are indicated 
in Table XI. The optical trigger which 
indicates blade position is located at 
X = 3.284 cms, r = 4.556 cms, and 9 = 52.9®. 


TABLE XI - 

■ IMPELLER 

AND 

INLET TRANSDUCER 

LOCATIONS 


Axial 


Radial 

Angular 

Number 

Distance, 

X 

Distance, r 

Location , 6 


(cms ) 


(cms) 

(degrees) 

K1 

-0.64 


4.40 

88.4 

K2 

0 


4.40 

67.2 

K3 

0.51 


4.40 

78.9 

K4 

1.02 


4.40 

92.5 

K5 

1.65 


4.40 

64.9 

K6 

2.18 


4.41 

80.6 

K7 

3.28 


4.56 

74.1 

K8 

5.24 


5.65 

87.0 

K9 

5.66 


6.31 

63.6 

KIO 

5.66 


6.31 

74.2 




X = 3.28 cms 


Optical 

trigger located at r = 4.56 cms 





9 = 52.9° 



b) Diffuser Locations - See Figure 38. 

There are eight (8) semi-conductor transducers 
located in the front cover and spanning the 
diffuser entrance and vane tip region. Five 
(5) of the transducers lie approximately along 
the central streamline with another 1-1/2 
throat width downstream of the diffuser throat. 
The remaining two (2) transducers are located 
near the vane tip. These taps are staggered 
circumferentially and located in several 
diffuser entry regions for clearance and so as 
not to interfere with the time-averaged static 
pressure probes. They are shown in one channel 
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in Figure 38 in order to show the region 
covered. The actual positions are given 
in Table XII. 


TABLE XII 

- DIFFUSER INLET TRANSDUCER LOCATIONS 

Tap Number 

Radial Distance, r 

Angular Location, 9 


(cms) 

(degrees) 

Kll 

8.45 

324.30 

K12 

8.61 

284.50 

K13 

8.84 

309.20 

K14 

8.99 

291.15 

K15 

9.18 

315,60 

K16 

8.51 

336.00 

K17 

8.46 

276.15 

K18 

9.98 

324.30 


Time-Averaged Total Pressure 

a) Inlet Plenum 

Two {2) wall static taps are used to obtain 
the plenum total pressure. 

b) Collector 

Five (5) wall static taps are spaced around 
the dump region of the collector. All of 
these taps should read the same pressure. A 
typical design of a static pressure tap is 
shown in Figure 40. 

c) Channel Diffuser Throat Taps 

Three (3) small total pressure probes are 
located in the front cover and approximately 
one throat width downstream of their respec- 
tive throats. This is well downstream of 
the position (0.468 x throat width) at which 
the channel area minus probe blockage equals 
the throat area. Figure 41 shows the probe 
design . 
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Figure 40~-Typic 0 1 Design of St a ti c Pressure 
Tap and Connector 
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Temperature Measurements 

a) Total (static) Temperature Inlet Plenum. 

A setup of three (3) Fe-Con thermocouples is 
used to make these measurements. 

b) Total (static) Temperature Collector. 

Five (5) Fe-Con shielded thermocouples are 
located around the colliictor. All of these 
thermocouples should read the same. 

Inlet Flow Surve y 

At the inlet to the inducer, a three-hole cobra 
probe is used to determine flow direction and 
static/total pressure. The cobra probe has a 
head width of 0.152 cm and is traversed manually 
in a radial line. Two axial locations at 0.625 cms 
and 2.54 cms ahead of the inducer lead edge will 
be used, with circumferential locations 345* and 
300®, respectively. 

Mechanical Instrumentation 


The steady-state and fluctuating components of the 
impeller axial motion at the tip region will be measured. 
This measurement is useful for diagnosing the mechanical 
behavior of the impeller, and for monitoring running 
clearances during system transients. 

A proximity type sensor has been found to be useful 
for this purpose. The sensor operates on the inductive 
proximity principle and therefore need not make electrical 
or mechanical contact with the rotating hardware. Mea- 
suring systems are available with operating ranges from 
about five (5) millimeters to about half a centimeter, 
and with a gain (output) or about 200 rav per 0.0025 cms. 

Because the probe is sensitive to any electrically 
conductive material in its vicinity, it should be cali- 
brated (output vs. displacement) with the same hardware 
in which it will ultimately be used. If this is not 
feasible, the calibration may be performed in a test fix- 
ture which replicates the materials and geometry in the 
area surrounding the probe at the anticipated installa- 
tion point. 
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This proximity sensing system is also sensitive to 
the temperature of the probe - a change of 11.30®C being 
approximately equivalent in output to a displacement 
change of 0.0025 cms. Thus, the probe should also be 
ca ■ ibrated over the expected temperature range. During 
actual use, a thermocouple or thermistor placed adjacent 
to the probe is used for accurate temperature indication. 
(See Figure 42.) 

Readout instrumentation for this instrument consists 
of a digital DC voltmeter for the steady state displace- 
ment and an oscilloscope for monitoring the fluctuating 
or transient behavior. 

The proximity measuring system consists of the 
probe, its signal conditioner and an 18 volt DC power 
supply. 


Readout Instrumentation 


All time-averaged static pressures will be read on 
Creare's pressure transducer system which includes pneu- 
matic switches, diaphragm-resistance bridge transducers 
with power supplies and a digital millivoltmeter. 

The dynamic component of the static pressures will 
be recorded with oscilloscope and oscillograph traces, 
depending on the frequency of the signal of interest. 

In addition to the dynamic pressure traces, a simul- 
taneous recording of the impeller position will be made. 
The impeller position measurement can be made using an 
optical pickup placed to look at the inducer blades. 

This will require prior calibration to correlate the blade 
position relative to the pickup with the output signal. 

Temperature measurements will be read out on a digi- 
tal millivoltmeter through an electrical switching system. 

The flow direction survey will use a water manometer 
for nulling the cobra probe. Total pressure will be mea- 
sured with the time-averaging transducers. Flow angle 
and radial position are measured with conventional appara- 
tus available in-house. 
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MECHANICAL DESIGN 


Mechanical, material and stress analyses were 
carried out for the several compressor components in 
order to assure adequate operation of the compressor 
stage during low speed model testing in Creare's facil- 
ities and to meet contract specifications for possible 
future testing of the compressor hardware at NASA’s 
facilities. 

Creare will conduct all testing on this compressor 
using the Low Speed Modeling facility located in Creare's 
laboratories. Thus, although the compressor aerodynamic 
and mechanical design has been carried out for 100% N^ 

(100% of design speed) = 75,000 RPM, testing of this com- 
pressor at equivalent 100% N^^ will occur at an actual 

rotational speed of approximately 52,500 RPM. While the 
mechanical rotational stress analysis has been conducted 
for speeds up to 120% of 75,000 RPM (120% of design speed 
in air) , the blade vibrational analyses and shaft criti- 
cal speed analyses were carried out for the lower LSM 
operating speeds. 

Basic LSM Rig Description 


A diagramatic view of the LSM rig showing the prin- 
cipal loop components is given in Figure 43. 

The compressor impeller is designed to mount on 
Creare's high speed gear box. The high speed rotating 
assembly is of overhung design. The entire compressor 
assembly is supported from a massive inlet plenum; the 
basic support structure of the assembly is a large back 
plate which carries the gear box and from which the dif- 
fuser and cover components are supported. The inlet is 
compliantly coupled to the plenum. 

This configuration has been used in past compressor 
testing at Creare with no mechanical difficulties encoun- 
tered. The entire non-rotating assembly, diffuser, 
impeller back plate and covers, has been designed to 
maintain critical dimensional tolerances, especially 
diffuser dimensions and impeller-to-cover running clear- 
ances at all speeds. This has been proved in compressor 
testing. 
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Figure 43 - Low Speed Modeling Test Facility 




The components of the compressor assembly are: 


~ high speed rotating components (shaft, 
bearings, impeller, and impeller mounting 
hardware) , 

- inlet, inlet centerbody and inlet-to-plenum 
hardware , 

“ impeller cover (shroud) , 

~ impeller back plate and seals, 

- diffuser components (diffuser front plate, 
diffuser back plate, diffuser vanes) , and 

- miscellaneous hardware (bolts, nuts, 0-rings, 
etc. ) , 

The major components of the compressor unit are 
displayed (in a meridional view layout) in Figure 44 
The components designed and fabricated on this program 
are indicated by the cross hatched parts of the figure. 

High Speed Rotating Components 

The principal analyses performed on the rotating 
hardware were : 

impeller material selection, 

impeller stress and vibration analyses, and 

impeller/shaft/bearing critical speed analyses. 

Material selection - Annealed titanium alloy 
(Ti-6A1~4V) forging was selected as the material for the 
impeller- The primary consideration in selecting titani- 
um was the maintanence of high strength/weight ratio at 
elevated temperatures. While several heat treated steels 
have superior absolute strength characteristics and sev- 
eral aluminum alloys are preferable on a low temperature 
strength/weight basis, titanium alloys are preferable for 
high strength application between 204 °C and 53 8 "C. 

The applicable material specification for selecting 
Ti-6A1-4V material is S.A.E. Specification AMS4928G. 

This specification defines minimum strength and ductility 
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(elongation and reduction in area) and exceeds the 
requireinents of the comparable military specification 
MIL-T-9047E. 

Impeller stress - The blade stresses (blade root) 
were calculated using simplified approximate calculation 
methods intended to be basically conservative and to 
verify that the final aerodynamic design would pass 
detailed mechanical analyses. Blade vibration modes and 
frequencies and disc stresses were calculated by the 
Boeing Company using NASTRAN and WASP finite element 
programs . 

Creare's BLSTRS program calculates centrifugal 
forces arising from the rotating blade mass. The root 

blade force F“/ ^ rto ^ pd (Vol/dr ) is computed over an incre- 

^t 

mental axial length. dVol/dr depends upon blade taper 
ratio and blade section profile. Output is force per 
unit axial length exerted by the blade on the disc. The 
program also uses blade root area to obtain the blade 
root stress. In the essentially radial portion of the 
impeller (i.e., near the exit), a subroutine resolves the 
centrifugal force into components parallel and normal to 
the backward leaning blade and uses the normal component 
to compute the root bending stress. Bending stress is 
based upon a simple cantilevered beam idealization. 

Figure 45 shows the blade root tensile, bending, and 
combined stresses for 90,000 RPM (120% design air speed) 
operation with the design blade contour. Table XIII 
displays the calculated blade root stresses. 

Wheel mass, moment of inertia, and estimated burst 
speed were calculated by Boeing. These values are dis- 
played in Table XIV. Burst speed was calculated at 
110,740 RPM for the disc at room temperature. A 107,000 
RPM (142.7% air) burst speed at an average disc tem- 
perature of 149°C was calculated. 

Blade vibration frequencies v/ere calculated via 
NASTRAN finite element program. The resulting Campbell 
diagram is shown in Figure 46. All blade frequencies 
within the anticipated operating speed range are inducer 
blade modes. The impeller tip modes are so high as to 
not be excited in the operating range. 
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TABLE XIII 

- IMPELLER BLADE 

ROOT STRESSES 

(90,000 

RPM) 

Calculation 

Station* 

Radial Stress 

,Newt , 0 > 

> 

Bending Stress 
( — 5— X 10 

lU 

MH**Stress 

X 10 » 

1 

2.34 

1.35 



2 

2.83 

0.87 



3 

3.30 

0.36 

2. 

.87 

4 

3.54 

0.94 

3. 

.40 

5 

3.92 

1.22 

4. 

, 03 

6 

4.05 

1.60 

4. 

,56 

7 

3.22 

1.84 

5. 

.09 

8 

2.68 

2.42 

5. 

.65 

9 

2.11 




10 

1.96 




11 

1.84 




12 

1.72 




13 

1.65 




14 

1.54 




15 

1.50 




16 

1.48 




17 

1.48 




18 

1.50 




19 

1,54 




20 

1.60 




* See Figure 45. 

** Mises-Hencky Stress. 
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TABLE XIV - PHYSICAL CHARACTERISTICS OF IMPELLER 

Polar Moment 




Weight 

(kg) 

of Inertia,! 
(kg. m^) 

Blade 


0.0118 

2.60 X 10”® 

Total 

for 19 Blades 

0.2240 

4.95 X lO"** 

Disc 


1.2250 

1.53 X 10-® 

Total 

Impeller 

1.4490 

2.02 X 10”® 


Minimum burst speed at 149 °C metal temperature is 107,000 
RPM. 

Maximum allowable speed for prolonged operation at 149 °C 
metal temperature is 80,000 RPM. 
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Figure 45— Def i n I tion of ImpeHer Stress Colculot ion Stations 
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Speed N, rpm (x 10“^) 


Figure 46— Impeller Campbell- Diagram 
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The Boeing Program WASP calculates disc stress, 
radial deflections, and axial deflections of the impeller. 
The principal results are listed in Table XV. Stresses 
and deflections scale as the square of the speed. Thus 
the peak centrifugal disc bore stress is 5.26 x 10® 
Newton/m^ at air design speed (75,000 RPM) and 2.56 x 10® 
Newton/m^ at LSM design speed (52,300 RPM). Impeller tip 
deflections at impeller design speed are 8.64 x 10“® m 
rearward and 1.397 x lO”** m radially. Blade stiffening 
of the disc, which may be significant in reducing axial 
deflections, is not handled by the Boeing WASP analysis. 

The maximum speed for long term operation with elas- 
tic (no creep) behavior is 80,000 RPM (107.7% air design 
speed, 153% LSM design speed) . A number of cautionary 
points must be noted with respect to disc stress and 
deflection. Bore diametric growth from centrifugal 
effects is 1.016 x 10"“* m at design speed in air 
(75,000 RPM) and 1.524 x lO"** m at 120% air design speed. 
Therefore, considerable care must be taken in mounting 
this impeller on its shaft if this design is to be spun 
in air at 75,000 RPM and over. Also, as implied by the 
7.58 X 10® Newton/m peak bore stress at 90,000 RPM, the 
bore will yield locally and stress relieve at 90,000 RPM 
(120% air design speed) . A residual compressive stress 
will then remain in the bore at rest. 120% air design 
speed operation should certainly be initially attempted 
only in a suitable spin pit facility with the bore to 
shaft fit suitably remade after spin. If significant 
long term service is contemplated for this wheel beyond 
the Creare experimental work, then low cycle fatigue 
analyses should also be performed. 

Critical speed analysis - The high speed rotating 
component assembly is shown in Figure 47. The critical 
speed analysis of this "overhung" assembly showed that 
the shaft first bending mode critical speed occurs at 
approximately 12,000 RPM. Two more critical speeds occur 
at 29,500 RPM and 30,700 RPM and these are the bearing 
bounce and rock modes, respectively. A critical also 
occurs at 33,500 RPM and is a system critical and not the 
shaft second bending mode. No other criticals were 
encountered up to 90,000 RPM, 

The first predicted critical occurs just below the 
speed where the first critical has been experimentally 
observed using a compressor rotating assembly with about 
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TABLE XV “ MAXIMUM DISC STRESS AND RADIAL AND 
AXIAL TIP DEFLECTIONS 


90,000 RPM 

(120% Bp 


75,000 RPM 52,000 RPM 


^^D air^ 


^^D LSM^ 


Maximum Disc Stress 

( Newto n io~®) 7.580 

Radial Impeller Tip 
Deflection (mm) 0.203 

Axial Impeller Tip 
Deflection 

(Rearward - (mm)) 0.122 


5.270 2.530 
0,142 0.069 
0.084 0.041 
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the same rotating parameters as this present design. The 
"swell” in vibration levels predicted between 29,500 and 
33,500 RPM have also been observed in past tests. 

The first critical represents no problem as the rig 
is run quickly through this critical speed at start-up 
and shut-down. The "swells" at 29,500 to 33,500 RPM also 
pose no problems since past experience has shown that 
these are relatively low amplitude, rotational vibrations 
and the rig can be run through these speeds (and even 
held at these speeds for some time) without appreciable 
impeller deflection. 


Inlet Hardware 


The inlet hardware is of standard design which has 
been used on previous testing of compressor hardware of 
similar pressure ratio and flow rate on the LSM rig. 

Since the inlet hardware is not subjected to severe 
thermal gradients during test, the only substantial 
causes of deflections and stresses in the inlet hardware 
are pressure differntials between the impeller, plenum 
and ambient. In the inlet region, the magnitude of 
forces from these pressure differentials are small and 
the design has proceeded by making the inlet hardware 
axisymmetric and sufficiently robust (thick sections of 
aluminum) to withstand these small forces with very small 
deflection. 


Impeller Cover 


The primary cause of impeller cover deflections are 
pressure loadings resulting from pressure differentials 
between inside and ambient pressures. The calculation of 
the loading caused by aerodynamic forces during rig oper- 
ation at design conditions shows that the cover may be 
expected to deflect (move) axially forward less than 
2.54 X 10“® m. During LSM rig evacuation, the axial de- 
flection will be less than 1.016 x 10“^ m. These deflec- 
tions are inconsequential to the maintenance of compressor 
unit critical tolerances during rig operation. Stress 
levels created by these aerodynamic loadings are very 
small. The cover hardware is made from aluminum. 

Diffuser Components 

The diffuser assembly is a "sandwich" design wherein 
the diffuser vanes are clamped in place and mechanically 
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sealed between the diffuser front and back covers. The 
diffuser vane and diffuser front and back clamping 
(mating) surfaces are ground and polished to provide the 
proper clamping force load distribution over the vane 
surfaces and to provide mechanical sealing between vane 
suction and pressure sides. The diffuser clamping bolts 
pass through clearance holes provided in the diffuser 
vanes. The diffuser vanes are positioned in their proper 
locations using aligning pins. Stress levels imposed 
upon the diffuser components are extremely small. This 
includes stresses created by pressure differentials and 
differences in the thermal expansion coefficients of 
materials used in the diffuser. All diffuser components 
are made from aluminum except the diffuser vanes (also 
bolts and aligning pins) which are steel. 

The front diffuser plate has a cut-out section at 
the impeller exit radius for the LDV viewing window. The 
LDV window is made from quartz, ground and polished to 
provide a high quality optical surface for the LDV opti- 
cal studies. This window spans several diffuser vane 
tips and mechanically supports and aerodynamically seals 
across the diffuser vanes. Analyses and experimental 
bench tests have been made of this window configuration 
to ensure that any problems associated with thermal 
mechanical stresses, differential thermal expansion and 
resultant leakage between window quartz and diffuser 
cover metal, and mechanical clamping of the diffuser 
vanes with the quartz window can be successfully handled 
in compressor test operation. Check-out testing of the 
compressor hardware under design speed conditions on a 
companion program has demonstrated the efficacy of this 
window design. 
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HARDWARE FABRICATION AND INSPECTION 


All of the compressor hardware was subcontracted to 
and manufactured by outside vendors. 

Impeller 

Two impellers were manufactured by pantograph 
machining techniques. Inspection of the impeller hard- 
ware compared to Creare supplied specifications was con- 
ducted by the vendor. A check inspection of some of the 
critical impeller dimensions (impeller inducer radius, 
impeller exit radius, blade thickness, etc.) was carried 
out by Creare. Figure 48 shows the finished impeller. 

Impeller Cover 


The impeller cover was manufactured to the specified 
inside meridional contour and inspected by the vendor. 
Overall dimensional specifications were checked by Creare. 
Figure 49 shows the finished impeller cover. The two 
photographs display (i) the LDV inducer window port and 
(ii) the cover static pressure tap locations and unsteady, 
dynamic pressure transducer location holes. 

Diffuser 


Figure 50 shows the diffuser front and rear cover 
plates; each is shown with several of the diffuser vanes 
in place. The aluminum covers were machined with stan- 
dard machining techniques and inspected by Creare. The 
diffuser vanes (17 to a set) were machined using tape con- 
trolled techniques to ensure accurate replication of the 
vane dimensions, particularly the vane suction surface. 

The vanes were inspected by Creare. Visible on the front 
cover photograph is the LDV window opening. In the vane- 
less/semi-vaneless region of the cover are the static 
pressure taps and the several dynamic pressure transducer 
locations. Figure 51 shows a close-up photograph of a 
portion of the static pressure tap field in the vicinity 
of the vane lead edge. 

Dynamic Balance of Impeller Assembly 

The rotating impeller assembly (impeller, pre-stress 
hub, impeller shaft and bearings and tie bolt and locknut) 
was dynamically balanced. Final balance of this assembly 
was achieved to less than 1.037 x 10”® kg/m. 


118 




Figure 48 - Compressor Impeller 


Oii/GflVAlj PAajfl 
fOOfi 


119 













OEIGINAL PAGE IS 
OF POOR QUALITY 


Figure 51 
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Figure 52 shows the assembled compressor hardware 
prior to installation into the LSM test loop. 
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CONCLUSIONS 


A "state-of-the-art", high pressure ratio (pr =8), 
small mass flow {0.907 kg/s), centrifugal compressor stage 
has been designed, fabricated and mechanically checked out. 
This stage has been specifically configured to serve as a 
workaorse compressor for laser Doppler velocimeter (LDV) 
studies of flow' in the impeller inducer and diffuser (and 
cvencuallv in the impeller) . The requirement of making LDV 
neasurements in the compressor impacted the design in two 
ways; (1) it was desirable to maintain a constant inducer 
tic radius and (2) it was desirable to have a large blade 
spacing at exit in order to minimize the time to acquire 
data with the LDV and yet still sample particle velocities 
over a small percentage of the blade— to— blade spacing . The 
impeller has no splitters and thus a higher blade loading 
within the radial portion of the impeller than would other- 
wise be desired. An axial inducer configuration with con- 
stant inducer tip radius is used. The high blade loading 
analyzed for the impeller may result in a lower than pre- 
dicted efficiency although this is not anticipated. 

The predicted efficiency, total-to-static, at the 
design point is 771. A 12% choke- to- surge flow range is 
anticipated. The material and stress analyses show ade- 
quate strength characteristics for all components of the 
stage up to 120% of design speed in air (90,000 RPM) . 

The impeller hub bore stresses are expected to exceed the 
the material yield stress causing local yielding of the 
bore material at 120% design speed. Appropriate spin 
tesring of the impeller is therefore called for at this speed 
if the impeller is eventually operated in air. The present 
use of the impeller is under low speed modeling (LSM) condi- 
tions where design speed will be approximately 52,500 RPM. 

A mechanical checkout of the entire stage up to 
52,500 RPM in the LSM test loop has verified the mechani- 
cal critical speed analysis of the high speed, rotating 
impeller/shaft/bearing assembly. Checkout at design 
speed has also verified the efficacy of the LDV inducer and 
diffuser window mechanical designs. No problems have been 
encountered with sealing or thermal or mechanical stresses. 

The compressor hardware, including spare cover, impeller 
and diffuser hardware is presently available for the conduct 
of basic centrifugal compressor research studied of the 
rotating impeller inducer and impeller exit/diffuser entry 
aerodynamics using LDV technology developed under a 
companion program. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


125 



APPENDIX A 

CREARE CENTRIFUGAL COMPRESSOR MODEL NOMENCLATURE 
AND CREDAS INPUT/OUTPUT NOMENCLATURE 
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Figure 2, Station Convention and Stage Terminology. 
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WHERE a = DIFFUSER ANGLE - FROM RADIAL TO "CENTER 
LINE" DIRECTION 

0 - DIFFUSER VANE SPACING (ANGULAR) 

26 = CHANNEL DIVERGENCE ANGLE 

W = THROAT WIDTH 

L = DIFFUSER LENGTH (OF BOUNDED REGION) 

= IMPELLER TIP RADIUS 

- DIFFUSER VANE LEADING EDGE RADIUS 

= DIFFUSER THROAT RADIUS 

1 « LENGTH OP CONSTANT AREA SECTION AHEAD OF 

THROAT 

Figure 3 . Vane-Island Diffuser Geometry. 





NOMENCLATURE 


eff 


geom 

AR. 


AR 

e 

AS^ 

a 

a 

B 


equiv 


flow area (normal to mean velocity vector, specifi- 
cally defined) 

effective area = ’'one dimensional" cross-sectional 
flow area 

geometrical cross-sectional area 

impeller passage geometrical area ratio (inlet to 
outlet) 

equivalent impeller one dimensional flow area ratio 

channel diffuser throat aspect ratio (b/W for radial- 
plane divergence; W/b for meridional-plane divergence) 

acceleration 

speed of sound 

boundary layer blockage; 

A 


B = 1 - 


effective 


BF 


geometrical 

meridional depth of passage (normal to mean 
meridional velocity component) 

blockage factor: 


„„ _ effective 
~ A 

geometrical 


absolute velocity (relative to a Newtonian frame, 
e.g. compressor casing) 

wall friction coefficients 

T 


Gf = 


pW^/2g, 


where : 

W is measured relative to subject wall 
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i 


J 


pressure recovery coefficient; 

P “ Pref 
“ (Po-P>ref 

where ; 

measuring and reference states and stations must 
be specifically defined, 

specific heat at constant pressure 

torque coefficient, friction 

diffusion ratio; 


DR = 


energy 

impeller effectiveness; MR 2 /MR 2 £ 


force 

proportionality constant in Newton's Second Law 
F = Ma/g^ 

boundary layer shape factor = 6*/9 

boundary layer shape factor for incompressible flow 
before transformation to compressible flow 

static enthalpy/unit mass 

stagnation enthalpy /unit mass 


total enthalpy /unit mass (in a coordinate system 
rotating in Newtonian space in absence of electricity, 
magnetism, gravity and capillarity) ; 


^ ^ ^ - u^ 

^ 2g 0 
^o 


incidence angle of flow onto blades; 


(i = - 3) 

constant = 1; (778 ft-lb^/Btu) 
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k 


L 

LP 


S, 

M 

M 

m 


m 

MR2 


MR 


2i 


M, 


rel 


N 


N 


s 


P 

Pq 

pr 


ratio of specific heats 

diffuser centerline length (from throat to exit plane) 
impeller load parameter: 

LP = 

length of channel diffuser throat, see Figure 3 

Mach number 

mass 

distance along specified meridional trajectory (often 
streamline of potential solution) 

mass flow rate 

Mach number ratio = it/^rel 2 j 

~ ^rel It'^^^rel 2i 


2i 


= ideal impeller discharge relative Mach number for 
a given impeller geometry, m and N and for isen- 
tropic internal flow; discharge blockage B2 = 


only blade metal blockage, uniform discharge 


state and = 3 


b2 mean line 


shaft speed 
specific speed: 


N = 
s 




(Ah^) 

o 


374 


where ; 


N = Radians/sec: (rpm) 2 3 

= inlet flow = in m'^/sec; (ff^/sec) 

Ah = ideal stage enthalpy rise in Joules/Kg: 

static pressure 
stagnation pressure 
pressure ratio 

P*" “ P/Poo 
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Pip total pressure; a rotating coordinate property defined 
by (for perfect gas) : 

k 



q dynamic pressure; 



Q heat 

R radius ratio, 

R gas constant 

®C,®K degrees 
r radius 

r radial coordinate (see Figure 1) 

S distance between impeller blades: 

S = 27rr/Z 

S distance along streamline 

SP stability parameter; SP = (9pr/9m)/pr 

s entropy/unit mass 

T static temperature 

stagnation temperature 

total temperature (defined by and for a perfect gas) 

t blade thickness (in direction specified) 

tj^ blade (metal) thickness 

u impeller (metal) velocity 

W diffuser throat width (in radial plane) 

W relative velocity (in coordinate system rotating 
steadily in Newtonian space) 
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w 


X 




V 


X 

z 

a 

6 

Y 

6 * 

e 

E 

n 


0 

6 

20 

X 


^WL 


total shaft work per unit mass of fluid 
work input to impeller per unit mass of fluid 
slip velocity: 


^s ~ ^m2 ^b2 “ ^02 

axial coordinate (see Figure 1) 
number of blades 

absolute flow angle (see Figure 1) 

relative flow angle (see Figure 1) 

pitch angle in meridional plane (see Figure 1) 

boundary layer displacement thickness 

wake width 

effectiveness (for impeller = MR/MR^) 
efficiency: 


n = 


h ^ - h 
os oo 

W 


(measuring stations must be specifically defined) 

boundary layer momentum thickness 

tangential angular coordinate (see Figure 1) 

diffuser divergence angle 

swirl parameter: 

X = Co/C 
0^ r 

work input coefficient: 

W 


U 

Welliver Loss Coef ficient= (P^2*“^o4^ ^ ^^o2*~^2*^ 

vorticity or loss coefficient: 
fiGv = 


134 



irMi impeller-tip Mach number: 
irMi = 

p density 

stagnation density 
0 slip factor; 

0=1- Vg/U2 

T fluid shear stress 

T torque 

^ yaw angle (between velocity vector and null 

direction of flow angle measuring probe) 

impeller angular velocity 

Subscripts 


0 , 1 , 2 , 3 , — 
b 

bf 

CL 

c 

coll 

crit 

D 

dm 
ef f 
f 

geom 

h 

i 


stations in the stage (see Figure 2) 

blade property 

back flow 

centerline 

cover 

collector station (receiving volume 
after diffuser) 

critical value (must be specifically 
defined) 

diffuser 

(impeller) discharge mixing 

effective 

friction 

geometrical 

hub 

inlet or impeller 
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i 

IW 

IGV 

j 

le 

m 

o 

n 

P 

P 

r 

rd 

ref 

rel 

s 

s 

sep 

surf 

T 

t 

TS 

TT 

WL 

w 


ideal 

internal wake (loss) 
inlet guide vane property 
jet 

leading edge 
mean 

stagnation 
normal thickness 
pressure surface of blade 
polytropic 

radial component (see Figure 1) 

rear disc (friction) 

reference state or station (must 
be specifically defined) 

relative to impeller coordinates 

indicates that process follows an 
isentropic path 

suction surface 

flow separation value 

surface property 

total (see definitions of 
and p^) 

tip or throat 

total to static (efficiency) 
total to total (efficiency) 
Welliver loss 
wake 


136 



X 


y 

e 

Superscripts 

* 


•k^ 


I 


O. 


Miscellaneous 


upstream of shock (e.g. M } or axial 
component {see Figure 1) 

downstream of shock (e.g. M^) 

tangential component (see Figure 1) 


mixed out state (must be specifically 
defined) 

impeller discharge state including 
all work input from fluid angular 
momentum, cover and rear disc 
friction and back flow 

time-varying property 

instrument-indicated property 
mass-flow-averaged property 

time-averaged property 

area-averaged property 

vector quantity 


In 

log 

exp 

e 

f C) 

tan ^ { ) 

A 

d 

3 

TT 


natural logarithm 
base 10 logarithm 
exp a = e^ 

base of the natural logarithms (e = 2.71828...) 
function of ( ) 
inverse operator: 

a == tan mean tan a = A 
incremental (but finite) change 

incremental (but infinitesimal) change, total derivative 

incremental (but infinitesimal) change, partial 
derivative 

3.1416. . 
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